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Abstract 

In this paper, the twist-3 two parton corrections in charmless B — ► PP decays are shown to 
be factorizable under the QCD factorization approach. The factorizability of the twist-3 two 
parton corrections is constructed on the following findings. Under the energetic meson limit, the 
pseudoscalar distribution amplitude for a light pseudoscalar meson is allowed to be non-constant 
by the equations of motion for the quark. The non-constant pseudoscalar distribution amplitude 
is then used to regularize the end-point divergences in the hard spectator corrections at twist-3 
order. By retaining the momentum fraction variable of the spectator quark of the B meson in 
the propagators, the end-point divergence in the weak annihilation corrections at twist-3 order 
is resolved. The factorization of the 0(a s ) corrections under the two parton approximation is 
shown valid up-to 0(l/m.fe) . The hard scattering kernels of order 0(a s ) and 0{h.QCD/mb) are 
explicitly given and found to be infrared finite. The results are applied for making predictions for 
the branching ratios of B — > ttK decays. 
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I. INTRODUCTION 



The hadronic B decays are good places for testing our understanding of standard model 
(SM) and/or new physics (NP). Until the present time, the B factories have obtained a lot 
of remarkable results. For example, the sin (2/5) has been determined precisely from the 
measurements of the mixing-induced CP-asymmetry in the B — > J/ipK s decay [U EJ [31 H] , 
the branching ratios for B — > -kit, tiK decays have been measured with only few percent 
errors [3 EJ [7], the CP asymmetry for B — > ti + K~ has been measured precisely [SJ [9j [TO] , 
and the B — > 7r°7r° decay mode has been confirmed experimentally with unexpectedly large 
branching ratio 0, HU [321 EES]- In contrast, the theory suffers from large uncertainties from 
the nonperturbative dynamics of QCD, which are involved in the matrix elements contained 
in the decay amplitudes. If the theoretical uncertainties can be largely reduced by effective 
methods, then a deterministic demonstration on the CP mechanism, an optimal extraction 
of CKM parameters from the wealth of experimental data, and a clean separation of NP 
from SM could be derived from the investigations of hadronic B decays [H] . 

The major breakthrough in reducing the theoretical uncertainties comes from an observa- 
tion [151 EE] > which is based on a two loop analysis [17] , that the decay amplitudes for nonlep- 
tonic B decays can be factorized under the heavy quark mass infinity limit [151 US H3 HH US] • 
The factorization of the decay amplitudes means that the amplitudes can be separated into 
perturbatively calculable short distance functions and nonperturbatively incalculable long 
distance functions. The finding leads to the QCD factorization (QCDF) [T5], [T7 ] ITS]. Sim- 
ilar facts have also been found under other approaches, the soft collinear effective theory 



(SCET) [261 E3 EH EH EQl EU E21 [33] and the perturbative QCD (PQCD) factorization 
[3ll[35l[36l[37J[38l[39lIinillIlll21. 



The important ingredient of QCDF is the factorization theorem [43] . which has been 
one successful method in the studies of hard scattering processes, such as deep inelastic 
scattering, Drell-Yan processes, etc. Unlike traditional hard scattering processes involving 
only one single energy scale, the hadronic B decays involve multiple energy scales: the 
electroweak scale \iew ~ Mw, the hard scale \in ~ nib, the hard collinear scale fine ~ 



yrribAQCD, and the soft scale fis ~ A-qcd- This multi-scale characteristic feature of hadronic 
B decays requires more efforts in the applications of factorization theorem. 

The factorization theorem for B — > MiM 2 processes with M 12 being light mesons has 
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been shown up-to leading order (LO) in the heavy mass expansion (or in the l/m& expansion) 
and next-to-leading order (NLO) in the loop expansion (or in the a s expansion). The a s 
means the strong coupling constant and m& is the b quark mass. Throughout this paper, the 
final state mesons in B — > M1M2 processes are restricted to be pseudoscalar light mesons, 
7r and K. The extension of the results derived in this paper to other light mesons can be 
done similarly. Under QCDF, the matrix element of a four fermion effective operator Qi at 
0(a s ) and LO in expansion can be expressed as 

(M 1 M 2 \Q i \B) = F? hh {q 2 ) jf 1 duTl 3 (u)<p M2 (u) + {M x <-> M 2 ) 

+ ^ d£ ^ du ^ dvT^^MO^MA^MM (1) 

JO Jo Jo 

where Ff Ml (q 2 ) with j = +,0 are the transition form factors , are the short distance 
hard scattering functions and (t>B,4>M 1 ,4>M 2 denote the long distance light cone distribution 
amplitudes (LCDAs) for the external mesons, B,Mi,M2, respectively. The function 
contains contributions of hard scale and the T? function contains contributions of hard 
and hard-collinear scales. The contributions below the soft scale are attributed to the 
meson LCDAs. The contributions between any two energy scales are calculated by means 
of the renormalization group method. The prediction power of factorization theorem given 
in Eq. (Ilj comes from the universality of the LCDAs. Once the LCDAs are determined 
experimentally in some processes, they can be applied for making predictions for other 
processes. 

However, the leading order predictions of QCDF based on Eq. Q can not consistently 
accommodate with experimental data for many decay processes. For example, the theo- 
retical predictions are about one half of the experimental data for the branching ratios of 
B — > ttK decays [201 EI]- For understanding the experiments, high order corrections to 
the factorization equation Eq. Q would be considered [44] . The higher order corrections 
mean higher loop corrections in the a s expansion or higher twist corrections in the 1/mj, 
expansion. As the higher order corrections are considered within QCDF formalism, it is an 
important issue to investigate whether the factorization of higher order corrections is still 
valid. 

So far, only partial results of order 0(a 2 s ) for T I,n have been calculated (IHl HHJ HT] HHl 
I49| [50] 151] 1521 [53] . These calculations showed that the 0(a 2 ) corrections seem preserving 
the factorization. On the other hand, the factorization of higher twist corrections is still 
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unclear. It was found that some high twist hard spectator and weak annihilation corrections 
spoil the factorization. The twist-3 hard spectator corrections [20] contain the divergent 
term 



where <f> p (u) is a twist-3 two particle pseudoscalar LCDA for a pseudoscalar light meson. 
Because 4> p {u) = 1 in the chiral limit [BH EE], X H diverges at u = 0. Similar end-point 
divergences also happen in the weak annihilation corrections, in which the divergent term is 



where <pp(u) is a twist-2 LCDA for a pseudoscalar light meson. These end-point divergences 
spoil the QCD factorization at the order of 0(a s ) and 0(l/m&). Since the Xh term are 
related to chirally enhanced corrections, which have a numerically large factor r x ~ 0(1), 
the power corrections can be equally important as the radiative corrections in charmless 
hadronic B decays [20]. From this respect, the extension of QCDF to subleading twist order 
is of some urgency [TH I2H] • 

However, the systematic generalization of QCDF to subleading twist order is still not 
available. For this reason, we restrict ourselves in this paper to only investigate the physics 
related to the non-factorizability of the Xh,a and search for a possible resolution to the 
related end-point divergences. In QCDF, the factorization is based on the collinear factor- 
ization scheme in which the partons participating in the hard scattering kernels are assumed 
to carry only collinear momenta. The hard scattering kernels are calculated by means of a 
leading twist collinear expansion method and perturbative QCD. The models for the LCDAs 
are submitted to nonperturbative theories, such as the QCD sum rules. Since the Xh,a terms 
are related to twist-3 contributions for B — > M\M.2 decays, the calculation scheme for the 
leading twist hard scattering kernels needs to be generalized. We will employ the collinear 
expansion (CE) proposed by Yeh [57] to calculate the relevant hard scattering functions 
up-to twist-3. The CE calculation scheme is a generalization of the leading twist collinear 
expansion for hard scattering processes jSU ESI EHJ EHJ EQ] ED E2J [63] . Similar extensions 
of the leading twist collinear expansion scheme for only two parton contributions have been 
proposed by Beneke and Neubert (the BN scheme) [20] and by Du, Yang and Zhu (the DYZ 
scheme) [23], El]. 




(2) 




(3) 
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Not only the calculation scheme needs to be generalized, but also the twist-3 LCDAs for 
the light mesons are required to be defined consistently to derive a factorization theorem 
at twist-3 order for B — > MiM 2 decays. For the Xjj term and other twist-3 two parton 
contributions for B — > M\M 2 decays, two twist-3 two particle LCDAs <fi p (u) and </><t(m) are 
involved. For a pseudoscalar light meson, there are three twist-3 LCDAs , the pseudoscalar 
LCDA 4> p (u), the pseudotensor LCDA 0<t(m), and the three particle LCDA <ps{u,u'). These 
three twist-3 LCDAs are related to each other by equations of motion (EOM)j6U [66]. Due 
to its small normalization factor, the three particle LCDA <f) 3 (u,u') is usually neglected 
in literature. However, this may not be a good approximation for some processes. For 
example, in the penguin dominated B decays, the tree level contributions from the three 
particle LCDA could be as large as the tree level contributions from the other two particle 
LCDAs [57]. 

In the approximation of neglecting the three particle contributions, <fi p {u) and ^> a {u) are 
determined by the following equations 



u ld(j) a {u) 1 

2 (MU) ~ 6^r ] = 6 

2 iMu) + 6^r ] = 6 



P\ U ) + ) = cMu) , (4) 



where u — 1 — u. The solutions are: (f> a (u) = 6uu and 4> p (u) = 1. We denote the solutions 
as the chiral limits of the 4> p (u) and (^(w). 

As explained above, the substitution of <fi p {u) = 1 into the Xh term results in an end-point 
divergence at u — 0. This end-point divergence may be due to the failure of the collinear 
factorization scheme, or the incorrect use of a model for the 4> p (u). In this paper, we intend 
to assume that the collinear factorization scheme is still applicable for the Xh term and to 
study a consistent model for the <i> p {u). This is different from the common viewpoint for this 
divergent problem as taken in the literature. In this respect, we found that the energetic 
meson limit is an important condition for solving the divergent problem in the Xh,a term. 
The energetic meson limit is defined as the limit at which the light meson's momentum 
becomes energetic. The energetic momentum means that the momentum contains a large 
component which is much larger than the light meson's mass and other components of 
the momentum. One light meson carries an energetic momentum is defined as an energetic 
meson. The energetic meson limit is the condition under which the leading twist factorization 
theorem Eq. ([T| has been shown to be valid. For example, the 0(a s ) radiative corrections 



for the Feynman diagrams as depicted in Fig. [IJa)-(d) can only be shown infrared finite by 
requiring the emitted meson M 2 to carry an energetic momentum such that the partons inside 
the M 2 meson can have large collinear momenta. The collinear (infrared) divergences cancel 
out in pairs between these four diagrams under the condition that the external partons to 
the radiative loops are collinear to their parent meson M 2 [TTl [18]. Because no any evidence 
shows that the higher twist contributions require different kinematics, we argue that, similar 
to the leading twist contributions, the twist-3 contributions (including the Xh,a) should be 
derived by using the same energetic meson limit, too. This argument has been used to 
derive the Xh term in [20] ■ Therefore, if the energetic meson limits for <p p {u) and <fi a {u) 
are different from the chiral limits, then one can hope to find a resolution to the divergent 
problem in Xh- 

However, the energetic meson limits for <p p {u) and <fio{u) have not been studied in liter- 
ature. We need to derive them in this paper. The details for derivations of the energetic 
meson limits for (p p (u) and (p a (u) will be given in Section [ll| Here, we briefly describe how 
the energetic meson limits for (/> p (w) and CT ( , u) can be different from the chiral limits. If 
the pseudoscalar light meson is energetic, then the 4>a(u) at u = 0(1) is found to be of 
0(A/E). The <f)' a {u) = 0(1) (the derivative of <p a {u)) and the <p p {u) = 0(1) at u = 0(1). 
The E denotes the energy of the energetic meson and A is of O(Aqcd)- It implies that the 
4> a (u) and the (f)' a {u) are of different order under the energetic meson limit, and they should 
be defined as different LCDAs. To keep both sides of Eq. ^ of the same order, the 4> a (u) 
should be dropped out and Eq. Q is further reduced to contain only the 4> p (u) and the 
(j>' a (u). The factor 1/6 associated with ^(u) in the Eq. ^ is only for normalization, it is 



then instructive to redefine the pseudotensor LCDA as that cf) a (u) defined in Eq. (23). By 



using the (f) a (u), the appropriate EOM for 4> p (u) in the energetic meson limit becomes 

M u ) = 4>v( u ) > ( 5 ) 

where the (f) p (u) and <fia(u) are defined as the energetic meson limits of 4> p (u) and 4> a (u), 
respectively. 

To have a better understanding of the above fact, we can define A0 (T ('u) = (<t> a {u) — 4> a (u)) 



as the difference between 4> a (u) and 4> a (u) by referring to Eq. (|8j) and Eq. (23). The (f) a (u) 
in Eq. |4]) should be A0 ct (m), and A0 ct (m) is of 0(A/E) in comparison to 4>' a (u) at u = 0(1). 
and should be identified as a twist-4 quantity in the energetic meson limit E 3> A. 
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Another interpretation is as following. According to Eq. (|8j), it is better to expand the 
coordinate variable zp in the spin projector for the <fi a (u) into its collinear and transverse 
parts, and to assign them by corresponding LCDAs, denoted by 4>' a {u) and (/) ff (u), respec- 
tively. To have a more transparent notation, we define 4>\{u) = 4>' a (u) and 0^(w) = (f) a (u), 
which respect the collinear degrees of freedom and the transverse degrees of freedom of the 
light meson state |M), respectively. These two <p\(u) and (j)^(u) are equally important at 
the condition E ~ A, which is applicable for a soft light meson. However, <j>\{u) and 4>^{u) 
become of different order as we boost the reference frame along the collinear direction of the 
light meson's momentum such that E ^> A, which is applicable for an energetic light meson. 
In the E ^> A reference frame, 4>„{u) becomes suppressed by a factor A/E than This 
is consistent with the parton model picture that the collinear degrees of freedom dominate 
over the transverse ones. In summary, we arrive at two consistent explanations for the same 
thing by noting that 4>\(u) cx (j) a (u) and 4>a( u ) °t 

To clearly clarify the source for the divergence in the Xh term, we also need to examine the 
calculation method related to the CT (w). The method proposed by Beneke and Neubert [20] 
is to separate the spin projector for the 4>a{u) into the collinear and transverse parts (the BN 
scheme). The collinear part is transformed into a derivative over the momentum fraction u 
and the (collinear) derivative is then defined to be applied on the (f) a (u). The transverse part 
is transformed into a momentum derivative in the transverse direction, and the (transverse) 
derivative is then defined to be applied on the hard scattering kernel. Alternatively, one 
can also let the collinear derivative applied on the hard scattering kernel and the transverse 
derivative applied on the </> CT (u). These two approaches are equivalent mathematically, but 
may result in different physical results. It is known as the projection ambiguity. To solve 
this ambiguity, Du, Yang and Zhu (DYZ) [23l 121] proposed that the whole momentum 
derivative should be applied on the hard scattering kernel (the DYZ scheme). As mentioned 
previously, the 4>a( u ) an d the <j>^{u) are of different magnitudes in the energetic meson limit, 
and their associated hard scattering kernels should be calculated separately. Since the above 
two calculation schemes did not consider the difference between 4>\(u) and 0^(w), the results 
calculated by these two schemes require further examinations. In this paper, we propose to 
employ the collinear expansion method to re-calculate the twist-3 two parton contributions. 



The collinear expansion method will be described in details in Section III 



For the weak non-singlet annihilation corrections, the end-point divergent problem is 
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more severe because it exists even for leading twist-2 LCDAs. This can be seen from the Xa 
term in Eq. ([3]). To solve this problem, we propose to retain the momentum fraction variable 
of the spectator quark of the B meson in the denominators of the parton propagators. In 
literature, the momentum fraction variable of the spectator quark of the B meson is always 
neglected. This is because the distribution function for the B meson is highly asymmetric 
such that the momentum fraction variable of the spectator quark of the B meson is of 
order A/m&. The reason to retain the momentum fraction variable of the spectator quark 
of the B meson is as following. The divergence in Xa term mainly arises as the parton 
propagators of the B meson become on-shell. To regularize the divergence in Xa term, we 
let the propagator be slightly off-shell by adding a term of 0(A 2 /ml) into it. After applying 
this for calculation, there are two same factors in the numerator and the denominator of the 
propagator of the spectator particle, respectively. As a result, these two factors cancel. The 
Xa in Eq. ^ becomes 

X A - / <%<I> B (Q / — — Mu) ■ (6) 
Jo Jo [u — £)U 

It is obvious that the original end-point divergence is regularized by the momentum fraction £ 
carried by the spectator particle. The only price we need to pay is to retain the 4>b{0 without 
integrating it out. From the resultant expression, the spectator particle is interpreted to 
carry a collinear momentum, although the spectator particle's momentum is soft. The above 
argument is valid for (V — A)(V ± A) and —2(5* — P)(S + P) operators. Therefore, the 
divergent problem associated with the Xa term is resolved. Similar fact for the factorizability 
of the annihilation contributions has been observed in [67J EH] 

The organization of this paper is as follows. In Section [TT| the energetic meson limits 



for (f) p (u) and (j) a {u) and related physics will be studied in details. In Section III, a simple 
introduction to the CE expansion scheme will be given first. The CE scheme is then applied 
to analyze the next-to-leading order radiative corrections. The factorization for the next-to- 
leading order radiative corrections is shown to be valid at 0(1/ E) or 0(l/ra&). In Section 



IV, the CE scheme is compared to the BN and DYZ schemes for the contributions related 



to the 4> a (u). In Section [V] the decay amplitudes at twist-3 order and at a s order will be 
recalculated under the CE scheme. The explicit expressions for the amplitudes for B — > PP 
decays, in which the final state PP means pseudoscalar light mesons, will be given in this 
Section, too. The predictions for the branching ratios of B — > irK decays will be also present. 
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The last Section is devoted for conclusions. The calculation details for the twist-3 0(a s ) 
vertex and penguin corrections are given in Appendix A and B. 

II. THE ENERGETIC MESON LIMIT AND THE LIGHT CONE DISTRIBUTION 
AMPLITUDES FOR LIGHT PSEUDOSCALAR MESONS 

In the B — > M\M.2 decays with M\$ being light mesons, the momenta P\ and P2 of the 
Mi and M 2 mesons can have a component being much larger than the other components 
and the meson mass. In the rest mass frame of the decaying B meson, the momentum 
conservation P B = Pi + P2 and the smallness of the meson masses m Ml ,m M2 <C M B lead 
to Pi ■ P 2 = (M 2 - m\ h - m 2 M2 )/2 ~ M 2 B /2 and Pf = (P+, P 2 " P 1± ) = (M B /y/2, 0, 0) and 
P2 = {P} i P2 t P21) — (0, M B /y/2, 0). The light mesons in this situation are identified as 
energetic mesons. The energetic mesons are not limited in the charmless hadronic B decays. 
The light mesons in the hadronic B — > DM decays with heavy-light final state mesons, 
or the exclusive hard scattering processes,e.g., 777* — > 7, 777* — > n, can all be considered 
as energetic mesons. For convenience, we parameterize the momentum P M of an energetic 
meson as 

p» = En" + ^|n^ , (7) 

in which E = M B j\pl and two light-like vectors n" = (n + ,n~ ,n±) = (1,0, 0j_) and n" = 
(n + ,n~ ,n±) = (0,1, 0j_) are introduced. The vectors n" and n" satisfy n 2 = n 2 = and 
n ■ n = 1. For later discussions, we define the chiral limit as E ~ rriM ~ ^-qcd, and the 
energetic meson limit asE> thm, ^qcd- 

The meson distribution amplitudes with a given twist order are defined as the momentum 
fraction distributions of partons in a particular Fock state of a meson. Up to twist-3, the 
distribution amplitudes of a pseudoscalar meson M are defined by the matrix element of 
nonlocal operators [64J 

(M(P) |g(0)[0;z]g(*)|0>|, 2=0 
= -i f -f jf 1 due^i^P.Mu) + /if 75 Up(u) - a^P^^y (8) 

where [0; z] denotes the Wilson line for preserving the gauge invariance of the matrix element. 
P is the meson momentum, fu the decay constant of the meson M, and /if = m 2 M /{m q + 
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mq) with tum, mq and rriq being the meson's mass and the quark and anti-quark current 
masses. The twist-2 distribution amplitude (f)p{u) and the twist-3 two particle distribution 
amplitudes <fi p (u) and <fi a {u) have been studied, in the chiral limit, based on nonlocal product 
expansion and conformal expansion. In addition, the equations of motion of on-shell quarks 
in the meson were used to obtain two differential-integral relations between the twist-3 two 
particle distribution amplitudes 4> p (u) and 4>„{u) and the twist-3 three particle distribution 
amplitude <f)^{u,u'). The differential-integral relations can be solved by means of moment. 
4> p {u) and 4> a {u) are determined by <p 3 (u,u') |64J. 

Some simplifications can be obtained by neglecting the twist-3 three particle distribution 
amplitude </>3 (u,u') due to its small normalization constant. The two differential-integral 
relations are then reduced to the relations as shown in Eq. Q. We identify the solutions 
to Eq. (jij) as the chiral limits 4> p (u) and 4> c a (u) of 4> p (u) and 4> a (u), in which 4> p (u) = 1 and 
(jf a {u) = 6uu with u = 1 — u. 

Referring to Eq. (pi), the Xh term from the hard spectator diagrams contains the partonic 
part 1/u and the hadronic part 4> p (u). As the chiral limit solution 4>p{u) is substituted into 
the hadronic part <f) p (u) of the Xh term, an end-point divergence arises as u — * 0. The 
end-point divergence spoils the factorization at the twist-3 order for B — > M\Mi processes. 
In literature, the common viewpoint is to identify that the breakdown of the factorization 
for the Xh term is due to the failure of the factorization scheme, the QCD factorization. 
In this paper, we propose to take another viewpoint that the source for the breakdown of 
the factorization for the Xh term may be due to the use of the chiral limit solution for the 
<f) p {u). Our consideration is the following. The partonic part of the Xh term is derived 
under the assumption that the external meson is taken in its energetic limit [20] . while the 
hadronic part is using the chiral limit solution for the relevant distribution amplitude. If 
the energetic limit solution for 4> p (u) can be different from the chiral limit solution, then we 
can expect to find a resolution for the end-point divergence. 

To see whether the above argument is correct, it is necessary to find out the energetic 
limits, 4> p (u) and 4> a {u), of the 4> p {u) and 0<j(m). Since the energetic limits have not been 
studied, we will derive them in this paper. The first main result of this paper is to show 
that the energetic limits, <p p (u) and 0<t(w), indeed exist and are different from the chiral 
limits. We will employ a simplified method to derive similar relations to those in Eq. Q for 
4> p {u) and 4> a (u). According to the relations for 4> p {u) and 0o-(m), a non-constant solution 



fO 



for the 4> p (u) will be obtained. We will also develop a expansion scheme consistent with 
assumption for taking the energetic limit for the light mesons. The expansion scheme is the 



second main result of this paper and will be given in Section |III| By using the expansion 
scheme to re-derive the Xh term and substituting the non-constant solution for the <j) p (u) 
into the hadronic part of the Xh term, the end-point divergence is then resolved. 

In the following derivation, only the asymptotic solutions for 4> p (u) and 4> a (u) are consid- 
ered. The complete solutions for the 4> p (u) and 4> a (u) by following the traditional approaches 
[OH E5J EH] will be given in another place. Usually, the twist-2 LCDA 4>p[u) is defined as 
the probability of the transition of the meson M into the q(u)q(l — u) pair at zero transverse 
distance. We assume that the same conditions are also applicable for the <j> p (u) and CT (u). 
To define the 4> p {u) and 4> a (u), we start from Eq. (Ji). For convenience, the coordinate z in 
the matrix element can be parameterized under the energetic limit as the following 



z 



z 2 ,E_ X 



Z±f3 



(9) 



2A '" M ' E 

where the variable A is for boost invariance in the collinear direction and z± is assumed to 
be of order 0(1/E). It is noted that the A is required to be large to insure that the collinear 
component z M = X/En fl dominates. This requirement for A is consistent with the condition 
for an energetic meson in the highly boost frame. The quark field q(z) in Eq. ^ can be 
expanded with respect to z^ = X/En fl as 

dq(z 



q(z) = q{z) + 



dz^ 



[z-zY + 



(10) 



where dots means terms of order with n > 2. By using Eqs. and Q, the spin 

projector [P a , zp] associated with (f) a (u) defined in Eq. d8]) can be written as 



[P a ,zp] = X[n a ,n p ] - 



m2 M Z ± 

4A 
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mi 



n a , np\ 



■- -| lit T\/T i- -i 

+E[n a , z ±l3 \ + 2^-ln a , zxp\ ■ 



(11) 



By substituting Eqs. (10) and (11), into Eq. (|8j), we arrive at the following identity up-to 
0(z ± ) 



(M{P)\q{0)a a ^ 5 q{X/En)\0) + (M(P)\q(0)a a ^ 5 



dq(z) 
dzf* 



|0> 



-if M ^ l-i[n a ,np] J^due^\ d -^ + E[n ai z^} f du^^- 
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It 



u=l 



+*[n Q ,n^ A ^j } (12) 

where the P ■ z in the phase factor exp(iuP ■ z) has been approximated to be A, and the 
terms proportional to Tom or 2;^ have been neglected. An integration by parts has been 



used to obtain the first and third terms in the right hand side of Eq. (12). According to 



the convention of [T8j [20], we let 4>' a (u) and 4> a (u) correspond to the collinear part and the 
transverse part of the spin projector, respectively. Comparing both sides of Eq. (12), we 
choose the following identities according to the order of the z± factor 

(M(P) \q(0)^q(\/En)\0) = -f M ^[n a ,n^ Jj^due^^M - (e^u))^} , 

(13) 

(M(P) |g(0)a Q/37 5^g(A/ J En)|0)< = -if M ^E[n a , z ±/3 ] [' due mx ^P~ ■ ( 14 ) 

A Jo 



From Eqs. (13), (14) and (ph, we observe the following issues needed for further exami- 



nations in the energetic limit: 

• The meaning of the momentum fraction u: according to the parton model, the mo- 
mentum fraction u is defined to be the ratio of the collinear parts of the parton's 
momentum and the meson's momentum. That is, if we let k denote the parton's 
momentum, then the momentum fraction is defined as u = k + /P + , where we have 
used the u to distinguish from the u. However, the momentum fraction u in Eq. ^ 
is not defined as the u. On the other hand, the u in Eq. ^ can be interpreted as 
the fraction of the whole momentum P. This can be seen from the exp(iuP ■ z) in 
Eq. (|8]). If the parton momentum is denoted as k, then u = k ■ z / P ■ z. By using the 
parameterizations for P and z, the u is expanded as 

u = — = ^ + k ■ n f (15) 

P ■ z \ 1 m li z ± v ; 

where the parton momentum fc M has been parameterized as 

k 2 

k^ = k-nn tl + -^n>* + k± . (16) 
2k -n K J 

The momentum fraction u carried by the collinear partons should be defined under 

the limit \z±\ — > 0, u = lim\ z± \_>o u(z±) = k ■ n/P ■ n. The fraction u in Eqs. (13) and 



(14) should be interpreted as it. In the following text, the u is always interpreted as u 



to simplify the notations. 
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The boundary condition: the boundary terms in Eq. ( 13 ) implies the following equation 



(0 CT (l)-e a CT (O) 



. 



(17) 



Although the solution 4> a (l) = 0o-(O) = can satisfy the above equation, but the end- 
point u — > 1 behavior of (f) a {u) (, or, by translation invariance, the end-point u — > 
behavior of <f><r(u)) is related to A. Since the A depends on the reference frame, this is 
in contradiction with the universality assumption for the </> CT (w). 



The definition for the <fi a (u): take a differentiation on both sides of Eq. (14) in zj_ as 



(M(P ) | v( )tr n r, A> < A/£'w M A/£'w ) 1 0) - -if.u(iME[r, n .<! ±lll } / due iaX< ^\ll 



We observe from Eqs. (13) and (18) that the same 4> a {u) corresponds to two different 



matrix elements in the energetic limit. There arise confusions which definition, Eq. ( 13 ) 



or Eq. (18) should be used for <j) a (u). 



The equations of motion: to avoid the above confusions in the definition for 

we suggest to re-define 4>%{u) = ^(u) and 4>^{u) = 4> a (u) according to Eq. (13) and 



Eq. (18), respectively. Because the boundary term in Eq. (13) may not vanish, the 



boundary term is defined to be absorbed by 4>\(u). By using (j)\{u) and 4>^(u), the 
first equation in Eq. ^ can be written as 



(19) 



To find out the energetic limit of the above equation, it is instructive to rewrite <p\{u) 



and as the following expressions by taking Fourier transformations for Eqs. (13) 



and (18) 



u 



u 



dX 



-iuX 



(M(P)|g(0)a^7 5 n Q n /3 g(A/En)|0) 



(20) 



dX e 



-iuX 



f M ^Jo 2tt E 



(M(P)\q(0)^7 5 n a d^X/En)q(X/En)\0) . (21) 



It is seen that, in Eq. (21), there is a large factor 1/E for 0^(w), which 
is of short distance. The transverse derivative d± in the matrix element 
(M(P)\q(0)a a /3'j5n a dj_(X/ En)q(X/ En)\0) corresponds to the transverse momentum kj_ 
for the quarks in the meson, which is of order 0(A). Therefore, 4>^(u) is of order 
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0(A/E) for u = 0(1) in the energetic limit. On the other hand, 4>\(u) and (j) p (u) are 
of order 0(1) for u = 0(1) in the energetic limit. The power counting for <p\(u) 



and 4> p (u) implies that Eq. (19) is reduced to 4> p (u) = 4>\(u)/Q in the energetic limit. 



Based on the above discussions, we arrive at the following definitions for the <fi p (u) and 
(u) 

(M(P)\q(0)i^ 5 q(X/En)\0) = f M $ £ due inX ^ p (u) , (22) 
(M(P)\q(Q)a a ^q(\n/E)\Q) = -f M ^[n a ,n p ] C due lUX ^(u) . (23) 



We now show that the equations of motion for <fi p (u) and (fi a (u) in the energetic meson 
limit is 

Mu) = (f> a (u) . (24) 
Let's start from the following equation 

poo (j \ p -iuX 

= / 7, — <M|g(0) 75 m ■ d(Xn/E)q(Xn/E)\0) - f Mt i% . (25) 

To arrive at the above equation, we have assumed the boundary conditions 

lim e- lUX (M\q(0)l5q(\n/E)\0) =0 , lim e-™ A (M|g(0) T5 g(An/£)|0) = -if M ^ , (26) 

X^oo X— >0 

where the first holds due to the large fluctuation in e~ luX under A — > oo, and the second is 
the normalization condition. By using the identity n ■ (id) = — (T a 0n a (d^) and the fact 
that 

roo AX p -iuX 

— — (M\q(0)w a pn a dP(\n/E)q(\n/E)\0) = -ifu^^u)) ~ ful*% , 
we then obtain 

POO (j\ Jj 

ifMli M (u<P p (u))= J Q — e-^ x (M\m^ym/En)q(Xn)\0) + zf M ^ I (u ( pM) , (27) 



where 4> a (u) has been defined in Eq. (23). 
For latter uses, we use the identity 



l5(Taf3 = ^e af3riX a vX (28) 
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to re-express 4> a {u) in a form as 

f M ^4>M = ~/ o |V™ A <M|g(0)ef a aP q(Xn/E)\0) , (29) 

in which e°f = e al3lX n 7 n\. This form for 4> a (u) is convenient for our following calculations 
of the hard scattering kernels. The equations of motion for the q in the energetic limit is 

i$(\n/E)q(\n/E) = + 0(1/E 2 ) , 

which is valid up-to 0(1 /E 2 ). We finally arrive at Eq. 

By using (j) a (u), there is no the so-called projection ambiguity in the calculations of the 
hard scattering kernels, because there are no coordinate variables in the spin projector of 
4>(t(u). In the energetic meson limit, the ambiguity that the same is defined for different 

components with different magnitudes of the spin projector has been avoided. We conclude 
that 4> p (u) and 4> a (u) are more appropriate than <p p (u) and <p c a (u) for uses in the calculation 
of the Xh term and similar twist-3 contributions. 

The model for 4> p (u) can be obtained by solving Eq. (fsj) . Since the partons in the two 
parton Fock state of a light pseudoscalar meson can share the meson's momentum equally, we 
can assume that the parameterization models for 4> p (u) and <p(r(u) can have the asymptotic 
form 6u(l — u), for simplicity. Because now the <j> p (u) is no longer a constant, if the 4> p (u) 
in Xh can be replaced by 4> p (u), then the end-point divergent problem can be resolved. 
However, this requires us to examine the contributions related to Xh term whether they 
are completely coming from the <fi p (u) under the energetic meson limit. It is necessary to 
develop an appropriate calculation scheme for the twist-3 contributions. 

III. THE COLLINEAR EXPANSION 

As mentioned previously, the spin projector for (j) a (u) contains a coordinate variable z. 
This leads to a projection ambiguity problem. For comparison, we will first review what 
this problem is. Next, we will introduce a collinear expansion method for calculations of the 
hard scattering kernels of the one loop corrections to the matrix element. 
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A. The spin projection ambiguity and infrared divergences 

As mentioned in previous sections, the spin projector [P a , Zp\ for 4> a {u) contains a coordi- 
nate variable z. There are two methods for performing the calculations of the contributions 
from the <p a (u). The first method proposed by Beneke and Neubert (BN)[20] is to separate 
the z in [P a , Zp\ into its collinear and transverse parts as 



where k is assumed to be the momentum carried by the quark. The d/du is defined to be 
applied on <fi a (u), and the d/dk±p is interpreted to act on the associated hard scattering 
kernel for <p a (u). However, it is also legal to let the whole momentum derivative d/dkp 
applied on the hard scattering kernel. The latter method was proposed by Du, Yang, and 
Zhu (DYZ) [23"l [25 • However, both methods cannot avoid the infrared divergence in the 
Xh term. 

B. A preface to collinear expansion method 

In this section, the CE calculation scheme proposed by Yeh in [57] will be used for 
calculations of 0(a s ) contributions. The CE scheme is only applicable to calculate the 
decay amplitudes in the energetic meson limit. Under the CE scheme, the calculated twist-3 
contributions are interpreted to be composed of collinear partons. A higher twist, which is 
composed of collinear partons, is usually called a dynamical power correction [56J. There 
exist other types of power corrections, such as the power corrections from soft gluon or 
renormalon contributions. We identify these as non-partonic power corrections. For these 
non-partonic power corrections, the CE may not be applicable. To include these non-partonic 
power corrections within QCDF requires further assumptions. For example, the soft gluon 
contributions are better determined by the QCD sum rules or the lattice QCD. In this 
work, we only investigate how the dynamical power corrections can be calculated by the CE 
method. 

The idea of CE was first made by Polizer [58] . The systematical method was latter devel- 
oped by Ellis, Furmanski and Petronzio (EFP) [HUES]. Using CE, Ellis el at showed that, for 
DIS processes, the twist-4 power corrections (corrections suppressed by Q~ 2 with Q the rel- 
evant hard scale in the processes) can be factorized into its short distance and long distance 




(30) 
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parts as the factorization of the leading twist contributions. However, a parton interpreta- 
tion for the twist-4 power corrections is lost. To recover the parton model picture for the 
twist-4 power corrections, Qiu [51] introduced a Feynman-diagram approach to re-formula 
the EFP's method. In this Feynman-diagram language, a parton model interpretation for 
the twist-4 power corrections becomes trivial. 



C. Preliminary 

The organization of this section is as following. We will first describe how CE can be 
applied for tree level diagrams. We then apply CE to calculate one loop corrections to the 
matrix elements under the two parton approximation. The factorization of the amplitudes 
of the one loop corrections will be shown up-to twist-3 order. In the following, we will use 
1/E instead of using l/m^ to represent the twist-3 order. To be specific, we shall consider 
the decay processes of a B meson into two pseudoscalar light mesons. The decay processes 
involve three restrictedly ordered energy scales: the W boson mass fiw ~ Mw scale, the 
factorization scale /ip ~ m^, and the characteristic energy scale of nonperturbative QCD 
Aqcd- With the help of operator product expansion (OPE), the relevant \AB\ = 1 effective 
Hamiltonian is given by 

H = %£ v q {C x {tiQ\{Li) + cMQ q M) 

V ^ q=u,c 
10 

+ CMQM) - ut(c 77 (/i)Q 77 (/i) 

k=3 

+C 8G (fi)Q SG (fi))] + H.C. , (31) 

where v q = V^V^for b — ► d transition) or v q = V q bV* t (for b — > s transition) and Cj(^t) are 
the Wilson coefficients which have been evaluated to next-to-leading order approximation 
by means of perturbative QCD and renormalization group [70j [TTJ [721 E3- The four quark 
operators Qi are given by 

Ql = (u a b a )v-A(qpUp) V - A ,Qt = (c a b a ) V -A(qpC(])v~A , 
Q2 = (Uab/3)v-A(Qf3Ua)v-A , Q2 = (c a b/3)v-A(q/3C a )v-A , 

Q3 = (q a b a )v-Aj2(vMv~A ,Qa = {qpb a ) V - A Q2(^ a qp) V - A , 

q> q' 

Q5 = (qab a )v-A^2{q'(3q'p)v+A ,Qe = {qi3b a )v-A^2{q' a q'p)v+A , 

q> q> 
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3 3 

3 
2' 



3 3 

Q9 = ^(q a b a )v-A^ e <l'(yW^V-A ,QlO = ^(qpb a )v-Aj2 e <l'(vWf3)v-A 



(32) 



and 



Q 77 = -^m 6 (g a or^(l + 75 )6 Q )F / 



""b\Ha^ V- 1 - 1 1i> u a ± uv 1 

8n 2 
9_ 
8tt 2 



Q 8G = -^ mb (q Q a^(l + l5 )t a a ^)G% , (33) 



where Q\ and Q 2 are the tree operators, Qz — Qq the QCD penguin operators, Qj — Q\q 
the electroweak penguin operators, and Q-j^ and Qs,g the magnetic and chromo-magnetic 
penguin operators. 

The contributions between the scales of /ip and nw are attributed to the Wilson coef- 
ficients, and the contributions between the scales of fip and A.qcd are then included into 
the matrix elements of the operators. By choosing an appropriate regularization method 
for the infrared singularities, the Wilson coefficients can be calculated to be independent of 
the external states and can be factorized from the matrix element (Mi(p)M 2 (q)\Qi\B(PB)) ■ 
For an energetic M 2 meson in the B — > M\M 2 processes, the color transparency leads 
to the naive factorization [HH] for the matrix element such that the matrix element 
{Mi(p)M 2 (q)\Qi\B(P B )) can be written as a product of a transition form factor and a decay 
constant in the following way 

(M 1 (p)M 2 (q)\Q l \B(P B )} « (M 1 (p)\Ji\B(P B ))(M 2 (q)\4\0) + {M x ^ M 2 ) . (34) 

The form factors are defined as 

(M 1 (p)\q(0)l"m\B) 
= F^(q 2 )(P^ + j?) + [F B - Ah (q 2 ) - F^tf)] ™* ~f M ^ . (35) 

The form factors coincide as q 2 = 0, F+~* 1 (0) = F ( f ->Ml (0). Due to the conservation of 
the currents J[ 2 , the scale invariance of the decay amplitudes is broken under the naive 
factorization. The radiative corrections to the matrix elements are necessary to compensate 
the scale dependence of the matrix element [151 EE3 US]- 
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The bottom meson momentum Pb is denned in a light-cone coordinate frame such that 
it can be written as P% = (p^ + g M ) with two light-like vectors = (g + , q~ , q]_) = En^ and 
pi* _ (p + ,p~,p % ± ) = En^ in the plus and minus directions, respectively. The Mi meson is 
defined to receive the spectator quark of the bottom meson as its quark component. The M 2 
meson is defined as the emitted meson from the weak interaction vertex. We associate the 
momentum p M for the M\ meson and the momentum q^ for the M 2 meson. The following 
metric and antisymmetric tensors are useful in our calculations 

w t' = & ~ > K' = $ - » 

ef = ^ ? ef = e^n a np . (36) 

These tensors satisfy the following identities 

q,w», = > = , d£ = d£ = 2 , ef = ef = 2 . (37) 

For a collinear loop parton of the M 2 meson, it is convenient to parametrize its momentum 
Z M into its components proportional to the meson momentum g M , the light-cone vector n u , 
and the transversal directions 

Z" = £ra" + + . (38) 

2E 

For convenience, we further define the collinear component, Z M , the on shell component, Z£, 
and the off shell component, lg, of the momentum Z M as 

Z> = n • Zn" , 
Z 2 

Z L - Z + 2n-/ n +Z ± , 

Z^ = ra" . (39) 

5 2n-Z K ' 

For the M 1 meson, the momenta of the collinear loop partons flowing through the Mi meson 
can be parametrized in a similar way. We parametrize the collinear momentum in terms of p^ 
and n fl . Because the mass effects from the mesons' masses and the light quarks' masses are 
very small in the decay processes considered in this work, we shall neglect them completely 
in this paper. In this approximation, we let rriMi — 0, i — 1,2, and q 2 = 0. And the loop 
partons are assumed massless. 
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According to Eq. (38), a parton propagator can be separated into its long distance part 



and short distance part. If we write the loop parton propagator as 

dH 



/''<//■:) = J J^y 4 e uiy - Z) [F L (1) + F s (l)} 

F L (y,z) + F s (y,z) , (40) 



where 



FM = J^r , F M = . (41) 

r — ie 2n ■ I — xe 

The Fl(1) propagator corresponds to the long distance part of the propagator, since 
F L (y,z) oc 8(y — z). The Fg(l) propagator represents the short distance part because 
Fg(y, z) oc S(y — z). We now describe one important property of the long distance propaga- 
tor Fl(1). As Fl{1) contacts with a ^n^ component of a vertex 7 M in the parton amplitudes, 
the ^rz M vertex will extract one short distance propagator Fs(l) and one interaction vertex 
i^f v from the hadron amplitude in the following ways 

The momentum factor (I — l) u is then absorbed by the hadron amplitude due to the Ward 
identity [5U|55l|56]. The above identity is obtained by a simple manipulation [56]. We first 



insert an identity 1 = [J] 2 ) /I 2 into the left hand side of Eq. (42) and expand the f as fa + 1/s 
as the following 



Since {1/ L ) 2 = = (fs) 2 , we then obtain 



(43) 



(44) 



P P n I 2 I 2 a 
where the first term Ijjjjs in the right hand side of the above equation leads to a vanishing 
result as it contacts with its left hand side iI/l/1 2 term. The only contribution can only come 
from the second term 1/sI/l- We further expand the J/l in the following way 

Due to (jj 2 = tfi 2 = 0, the remaining result becomes 
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(a) W 

FIG. 1: (a)The type I diagram for B — > M\M 2 decays, (b) The type II diagram for B — » M\M 2 
decays. The central blobs with symbols H I,n represent the parton scattering functions. The initial 
state B meson is represented by the circle with a symbol B in the left hand side of each diagram. 
The final state M\ meson is represented by the circle with a symbol M\ in the right hand side 
of each diagram. The final state M 2 meson is represented by the circle with a symbol M 2 in the 
upper part of each diagram. The lines with a arrow are the fermion partons. 



By substituting the above back into Eq. (43), Eq. (42) is derived by noting that 



Un«m)i{i-i) a = uu 



Using Eq. (42), one can systematically include the effects from the non-collinearity and 
off-shell-ness of collinear partons external to the hard scattering function of a Feynman 
diagram. This property of the long distance propagator plays an important role in our 
following analysis. 

Let's first describe how the collinear expansion method can be applied for tree level 
amplitudes to recover the NF amplitudes. After the tree level analysis, we will investigate 
how CE method can be applied for the one loop amplitudes and show that the factorization 
theorem Eq. ([!]) up-to 0(1/ E). We begin from the diagrams of Fig. [l]for the B — ► M\M 2 
process. The parton interactions in the diagrams of Fig. [I] (a) and (b) are assumed by an 
effective four quark operator Qi = (q a T ib)(q b T iq c ) . In order to discuss the collinear expansion, 
we propose to formally write the transition matrix element (MiM 2 \Qi\B) involved in the 
amplitude for the diagrams of Fig. [I] (a) and (b) as the following expression 

+ ■■■, (45) 
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where F?~^ Ml (q 2 ) represents the transition form factor denoted by the bottom part of the 
diagram Fig. [I] (a). The scattering functions H I,n contain the parton interactions. The 
H 1 function represents four parton interactions with multiple radiative gluons. The H n 
function represents six parton interactions with multiple radiative gluons. The dots denote 
higher order terms which could contain interactions between the three parton Fock state 
of M 2 and the other partons from the M\ or B mesons. The contributions from the three 
parton Fock state of M\ or M 2 are neglected completely in this paper. The dots' terms are 
not shown explicitly. 

The diagrams in Fig. [I] (a) and (b) represent the processes that the initial state B un- 
dertakes a transition into the final state M\ by means of the b — » q\q 2 q 3 decays accompa- 
nying multiple radiative interactions. The partonic radiative interactions associated with 
b — > q\q 2 qz decays are collected into H 1 ' 11 . The probability for the transition of B into bq s 
pair is denoted by 4>b{Ib)- The q 2 q 3 pair produced from the interaction center, i.e. H I,n , 
then combine to form the M 2 meson after a long distance travel away from the interaction 
center. The probability for the transition of the q 2 q 3 pair into the M 2 is represented by 
4> M2 {Im 2 )- I n t ne diagram of Fig. [I] (a), the spectator quark q s of B combines with the q\ 
produced from H 1 to form the M\. The transition from B into M\ is represented by F B ^ Ml . 
In the diagram of Fig. \U (b), the q s gets involved in H 11 and then combines with q\ to form 



M\. The probability for the transition of q s q\ into M% is represented by (j) > 



Mi, 



We first discuss the diagram in Fig. [T] (a), which is represented by the first term of 



Eq. (45). In the expression of the first term of Eq. (45), the scattering kernel H*AIm 2 ) and 
the meson amplitude 4> M2 (Im 2 ) are correlated by the loop parton momentum Im 2 , the color 
indices, and the spin indices. The loop parton momentum Im 2 is defined to flow from the 
antiquark line to the quark line of M%. The expression for (f> 2 (Im 2 ) is defined as 

<P M2 (Im 2 ) = I d*ze- lI ^ z (M 2 (q)\q 2 (0)q 3 (z)\0) , (46) 



where Im 2 = q — hi 2 and the color and spin indices are not shown explicitly. The Tr 
denote traces over the color and spin indices. To complete the factorization of the first 



term of Eq. (45) into the short distance and long distance parts, we need to disentangle the 
correlations in the integration over Im 2 , the color indices, and the spin indices. 

In order to derive a factorization theorem similar to Eq. ([IJ, we propose to employ the 
expansion scheme made in (57]. The H^{Im 2 ) is expanded in a s and Im 2 ±- Similarly, the 
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H- 1 is expanded in a s and /m,±, i = 1,2. Namely, we first expand H^{Im 2 ) m at s as 

= 4 ( ° ) +4 (1) ^M 2 ) + -- - , (47) 
H^^Mn Im 2 ) — H- 1 1 (Zmij ^m 2 ) + • • • (48) 

where H 1 ^ is of order 0(oP s ) and -f^ 1 ^ and Hj 1 ^ are of order 0(q4). The dots are of order 



0(a™) with n > 2. Each term in Eq. (47) is then expanded in Im-i± or Im 2 ±- The expansion 
of and H*j in or Im 2 ± are left to latter discussions for one loop corrections. 



D. Tree level expansion 

The H-j ^ is independent of Im 2 and H 1 ^ = Ti5ij5 ab . The color factor 5 a b is then absorbed 



by cf) M2 . By substituting H 1 ^ back into the first term of Eq. (|45|), we arrive at 



Al = Ff^tf) I ^Tr[H!}°U M2 (lM 2 )] ■ (49) 



(2tt) 

We further use the integration transformation 



P 1 /" 1 pOO rf \ 

1 = / - / Ma ■n/q-n)= dv ^e^-h^-n/E) ^ 



/0 JO jO 27T 

to rewrite Ai as 

A a = Ff-* M '(q 2 ) £ dvTriHiP^iv)} . (51) 
The distribution amplitude <p M2 (v) is given by 

poo rl\r, 

^(v)= / ^ 2 e- A2 M2 (A 2 ) (52) 

JO Z7T 

where 

M2 (A 2 ) = J d 4 ze-^ z J |^ e ^-^- A2 ^(M 2 (P 2 )|g 6 (0)g c (z)|0) . (53) 
We further use the fact that the integrations over Im 2 and z result in a delta function 

/ dh I d ^^ lM2{z - X2nlE) = J d A z5^(z - X 2 n/E) . (54) 
By this, we write <fi M2 (v) as 

poo A \ _ in 

M2 M = y o ^e-^(M 2 (g)|g 2 (0)g 3 (A 2 -)|0), (55) 
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where v — 1 — v. To factorize the remaining spin indices in the trace, we make use of the 
Fierz identity 

1 1 

SijSki = ^((Y)ij(^)ki + {l^lbhAlbl^ki + (75)ij(75)fc« + ^{v^lbhAo^lsjki) (56) 

to obtain 

A, ~ F™tf) £ £ dv\Tr[H^Y' k ]Tr[Y' k r\v)\ , (57) 

where T' k T' k = ( 7 ")(7m), (75X75), |(^7 5 )(^7 5 ) for A; = 1, • • • , 4. To project 

onto the specific spin state of M 2 , we introduce the following definitions 

Tr[Y<p M2 {v)} = , 

Tr[Yl 5 <P M2 (v)} = -if M2 En^(v) , 

Tr[-f 5 (f)M 2 (v)} = -if M2 ^M 2 4>p 2 (v) , 

Tr[a^^{v)} = -f M2 fi M2 [n^n^(v) , (58) 

where <pp 2 {v) is the twist-2 LCDA, ^ 2 {v) and 0f 2 (w) are twist-3 LCDAs, and the LCDAs of 
twist order higher than three have been omitted. The final result depends on Dirac structure 
of the effective operator Qi. The Dirac structure can be Yi <g> Yi = (7^(1 — 75)) (g) (7^(1 ±75)) 
for (V-A)(V±A) operators, and T i ®T i = -2((1 - 75)) <g> ((1 + 75)) for -2(S-P)(S + P) 
operators. By substituting the Dirac matrices for the operators and using the definitions for 
the spin state of M 2 , the amplitude becomes 

J 

A^ S ~ P ^ = if M2 F^ (g 2 )m|rf jf 1 dv%* (v) , (59) 

where r^ 2 = 2m 2 M2 / {fhi ) {fh q2 + fhq z )) with rh b ,m g2 ,m q3 being current quark masses. In the 
above example, we have present how the collinear expansion method can be applied for 
tree level diagrams. As a necessary condition, it can recover the result obtained by naive 
factorization (NF). 

E. Factorizable one loop diagrams 

The 0(a s ) scattering function H-j(l M2 ) are classified into factorizable and non- 
factorizable parts 

4 (1) (/m 2 ) = H!}^ F (l M2 ) + H^' NF (l M2 ) . (60) 
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FIG. 2: The one loop factorizable diagrams. The external meson states are not shown. 






(c) (d) 





(e) (f) (g) 



FIG. 3: (a)-(d) The vertex diagrams, (e)-(g) The penguin diagrams. 

The H^' F (Im 2 ) corresponds to the the diagram in Fig. [2] (d). However, there are factoriz- 
able one loop diagrams as depicted in Fig. [2] (a)-(c), which are attributed to the B — > Mi 
form factors F?~* Mx . The sum of the first two diagrams Fig. [2] (a)-(b) for the form factor 
are dominated by soft gluons. Under QCD factorization approach, these leading soft con- 
tributions are defined to be absorbed by the physical form factors Fj 8- * 1 [T7] . The third 
factorizable diagram is to re-normalize the V — A current associated with the b — > q\ process. 
The fourth factorizable diagram is to re-normalize the V — A current associated with the 
M2 meson and the equations of motion. Since the V — A current is conserved, the V — A 
current and the equations of motion receive no renormalization. As a result, the factorizable 
diagrams lead to finite contributions under QCDF approach. 



F. Non-factorizable one loop diagrams and collinear expansion 

The H-j^' NF '{Im 2 ) corresponds to the non-factorizable one loop diagrams in Fig.pKa)-(g). 
According to the types of the diagrams, we write H*j' NF (Im 2 ) as 



H-j lS> ' Nf \Im 2 ) — Vij(lM 2 ) + Pij{l>M 2 ) + Psg,ij(lM 2 ) , (61) 
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where Vy denotes the vertex corrections from the diagrams in Fig. [3] (a)-(d), Py denotes 
the penguin corrections from the diagrams in Fig. [3] (e)-(f), and Ps g ,ij denotes the magnetic 
dipole moment corrections from the diagrams in Fig. [3] (g). Each type of corrections will be 
analyzed in following each subsection, respectively. 

We take the vertex correction Vij(lM 2 ) as an example to show the main feature of the 
application of the CE method to the one loop corrections. In the following, we omit the 
subscript M 2 in l M2 to simplify the notation. We write the amplitude for the diagrams in 
Fig. § (a)-(d) as 

A v ~ F^V) / ^nvm M (l)} ■ (62) 

In the above expression, the color and spin indices and the subscription of the form factor 
have been omitted for simplification. The irrelevant factors associated with the amplitude 
are also omitted for simplicity. To separate the collinear limiting V(l) of V(l) from the 
others, we make a Taylor expansion for V(l) with respect to I as 

dV(l) 



V(l) = V(l) + 



dl a 



(1-IT + ---, (63) 



i=i 

where dots are higher derivative terms. The first derivative term is shown for our latter 
comparison with the BN and DYZ schemes and its related contribution will not be considered 
in this work. By using 

Ot Ct — OL 

W a> = 9a' ~ n n a' , 

we can write 

(/ - i) a = w a j a ' . 

l a ' is then absorbed into 4> M2 (l) by Ward identity. With substitution of the first two terms 



of Eq. ( 63 ) back into Eq. (|62j) , the result appears 

Av ~ / (^i^^^^Wl + / ^^lUv,v)w^ a '(l,l)] + ■■■ , (64) 
where we have used the low energy theorems 

V(v) = V(l) , (65) 
dV(l) 



V a (v,v) = 



dl a 



(66) 
i=l 
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m, i ; b 




FIG. 4: The diagrams for V« (Z, Z'). For each diagram, there are similar diagrams with different 
gluon attachment. 



and have denned 



^ a '(l,l) = I d 4 ze~ u - z (M 2 (q)\q(0)id a '(z)q(z)\0) . 



(67) 



The parameterization Z M = vEn^ 1 has been used. The v has the meaning of momentum 
fraction carried by the partons. The collinear momenta I of the partons are defined to be 
parallel to the M 2 's momentum q in the direction. The functions V(v) and V a (v,v) have 
been introduced to emphasize the fact that, under the collinear limiting, V(l) and V a (l,l) 



and <j) M2 (l) are only correlated by v. Similar integral transformations to Eq. (50) can be 



A v ~\ dvTr[V{v)<l> M2 {v)} + [ dvTr[V a {v,v)wZ,<f>f 2a '{v,v)} + 



used to rewrite Eq. (64) as 
where <fi M2 (v) and 



-iXv 



n 



n 



(M 2 (q)\q 2 (0)zd a (A-) ?3 (A-)|0) 



(68) 



(69) 



have been used. 

For gauge invariance, we need to consider the diagrams as depicted in Fig. |4| We write 
the amplitude for the diagrams in Fig. [4] as 

dH r dH' 



(q 2 



(2tt) 4 7 (2tt) 



r Tr[V a (UW a (M')] 



(70) 
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where the function V a (l, V) represents 

V a (l, V) = V«(J, V) + vy>(j, V) + V^(l, I') + V^(l, I') + V^(l, V) 



(71) 



Each term in the right hand side of Eq. ( 71 ) corresponds to the diagrams as depicted in 



Fig. [4] (a)-(e), respectively. For each diagram in Fig. |4j there are other similar diagrams with 
different gluons' attachments. They are not shown explicitly. The amplitude (j)^ 2a (l,l') is 
defined as 



cf>^ a (l,l') = Jd'yJ dhe-^'- J >ye-^{MM\m{-)9^ a {y)T a q{zm . 



(72) 



In the following, we only consider the contributions related to V^ b \l, I'), which results in con- 
tributions related to (M 2 |g(0)i-D M (z')g(z)|0). The other contributions, i 7^ b, are related 
to the (M 2 \q{0)G flu {z')q(z)\0). The twist-3 contributions related to (M 2 \q{0)G' lu {z')q{z)\0) 
are left to our another preparing paper. 

The collinear limiting part of V^(l, I') is derived by using CE 



vV(i,i')=vV(V>)+ y: 



(k - kf + 



(73) 



By substituting Eq. (73) into Eq. (70) and neglecting the other terms except of V^ b \l, /'), 
we arrive at 

dH r d A V 



A' v ~ F B - M ^q 2 ^ 



:Tr[V^\iJ')^ a (l,n]. 



(2vr) 4 7 (2tt) 4 

Further more, it is convenient to rewrite A a (z') in <j)^ 2a (l,l') as 

A a (z') = n a n ■ A(z') + w«,A a '(z') . 



(74) 



(75) 



By substituting the above expansion for A a (z') into Eq. (74), we arrive at 

dH f dH' 



iB->M U 2 



(2tt) 4 7 (2tt) 4 
x{TT[V^(i,i')n a n a ^ a '(l } l')] 

+Tr[^)(f,f)<,0f a ' (/,/')]• 



(76) 



For covariant gauge, the first term of Eq. (76) can be transferred in the gauge phase 



factor and the second term of Eq. (|76[) can be combined with the second term of Eq. (64) 



(See below further explanations.). For physical gauge, such as n ■ A = 0, it is automatically 
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vanishing. In this work, we choose to use the covariant gauge in our following analysis. 
Under the covariant gauge, there are other contributions from diagrams with more partonic 
gluons of M 2 . The collinear limitings of these contributions are equally important and should 
be considered. However, they can be shown to contribute to the gauge phase factor. In the 
following, we exhibit this fact by only considering one partonic gluon case, the first term of 



Eq. (76). We now explain how the first term of Eq. (76) can be transferred into a gauge 



phase factor. By using the identity 



vw(i,i')n a = vw(i,i') 



n ■ k 



k=(i'-i) 



and the Feynman identity 



k a vi b \i,i') = v{i)-v{i') 



we can rewrite the first term of Eq. ( 76 ) in the following form 

dH 



dH' 



(27r) 4 7 (2vr) 4 



(77) 



(78) 



(79) 



where 



<t>n-W) = Jd A yj d*ze lk <y- n ^e- llz (M 2 (q)\q(0)(-w) J q d V n ■ A a (y)T a q(z)\0) . (80) 
in which the identity 

dr]e- ivn - k (81) 



(82) 



n ■ k — ie Jo 

has been used. Completing the momentum and coordinate integrals, we obtain 
J dv'Tr[V^(v)^ A (v)} - J duTr[V^{v')^.\{v')\ , 



where 



00 d\ 



-ivX i 



(M 2 (q)\q(0)(-ig) / drjn ■ A a (r]n/E)T a q(Xn/E)\0) 



(83) 



/O 27T JO 

Since V a (l, 1) and V^Q,, V) have similar structures, this enable us to add up Ay and A' v 
to obtain 

Ay + Ay ~ F B ^ Ml (q 2 ) 



(27T) 



+ 



dH r dH' 



(2vr) 4 7 (2vr) 



(84) 
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where 

4>D 2a (l, = / d'yj d i ze-^- l ^e- llz {M 2 {q)\q{^D a {y)q{zm . (85) 
with iD a = id a — gA a . By using the integral transformations for I and /', we arrive at 

+ jT 1 cfo jf 1 d^Tr^ 6 ^, </)<^ 2Q >, «')]) , (86) 

where 

V(v) = V{1) , 
(v,v f ) ee V?\tJ), 

POO rl\ 

<j> M *(v) = / ~e-™ x (M 2 (q)\q(0)<l(*n/E)\0) , 

JO 27T 

^'M') = Hp- r ^e-^-^e-^\M 2 (q)\q(0)^ a '(w/E)q(\n/E)\0) . (87) 

JO 27T JO 

The contributions associated with the V^ b '(v,v') are at least of 0(1/ E 2 ) and will be ne- 
glected. In the above, we have written the term related to V£ > '(v,v'). It is given here for 
comparison of the CE scheme with the BN and DYZ schemes. We have introduced how the 
CE is applied to derive the contributions related to the collinear limiting part V(v) and how 
the CE can separate different contributions related to different number of collinear partons 
of Mi- Because we are only concerning the twist-3 contributions, our remaining task is to 
show that V(v) is infrared finite up-to 0(1/ E). This is given below. 



G. Collinear expansion for vertex corrections 

The relevant term in the amplitude for the four diagrams depicted in Fig. [3] (a)-(d) is 
written as 

d% 
WY 



A Q< „ F B->M 1{g2) I _2_ Tr[ ^. (/)0 M 2(/)] (g8) 



where 
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where we have employed equations of motion for the b and q\ quarks for the terms inside 
the first trace bracket. The virtual gluon's momentum is represented by k. We first perform 
the following expansion for Vij(l) to derive its collinear limiting part V(l) 

dV 



W) = V{1) + ^ 



(I - l) Q + • • • . (90) 



1=1 



The virtual gluons could be hard k^ ~ (E,E,E), soft kg ~ (A, A, A) , or collinear to 
as kg ~ (E,X 2 /E,X) or collinear to p^ as k c , ~ (X 2 /E,E,X). Therefore, we divide the k 
integral into three regions corresponding to the soft ks, the collinear k c or kc, and the hard 
k H . 

To analyze the infrared structure of V^(l), we define the soft and collinear limiting parts 
of Vij(l), in which k are set as soft ks or collinear kc or kc, as V^ S {1) and V^ c (t) and 
^\c'(0> respectively 



1. Soft finiteness 



We first write Tr[V^fl(l)(f) M2 (l)} in its explicit form as 



In order to find the 0(1/E) contributions in the above expression, we have written the full 
part (i.e. I + k s or I + k s ) of the internal parton propagators in the upper part of the 
diagrams in Fig. [3] (a)-(d). Take an example for explanation, we consider the propagator 
for the left internal parton propagator in the upper part of the first diagram in Fig. [3] (a). 
Because k s <C I ■ ks, the propagator appears as 

- * (92) 



(l + k s ) 2 21 -k s In- 1 2n ■ I 
where i'f is an unit vector in the transverse directions, t 2 , =1. Since the denominator of 



the first term in Eq. (92) is of O(EA), and the denominators of the last two terms are of 
0(E), we obtain the first term is of 0(1 /A) and the last two terms are of 0(1/ E). Since 
E 3> A, the leading contribution comes from the first term and the last two terms are power 
suppressed as 0(A/E) in comparison with the leading first term. Since the last two terms 
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are independent of ks, they can be decomposed from the loop momentum integrations over 
ks- 

After we have separated the leading and subleading terms, where the latter are power 
and soft suppressed as 0(A/E) than the former, we now show that the leading part of Vs 



gives a vanishing result. By observing Eq. (92), we can see that the relevant terms of Vs are 



^WT S ~ ^Ts m{ - r vk f ' ~ ' (93) 

The 7° in the second trace term can be 7° = fin a , fin a , or, d" j a . If 7" = fin a or fin a , the 
second trace vanishes as 

21 ■ ks 11 ' fc 5 

= Tr[(^-/-f 4 -f 4 ^L^)^(/)] 
21 -k s 21- k s 

= , (94) 

where the first line vanishes due to fil oc n 2 = 0, and the second line vanishes since fil/ = vE, 
fif — vE and 



2n-l ' 2n-l 
= -(^-^ = 0. (95) 

If 7 a = d" 7°* , the contraction of d^/ with the first trace gives a vanishing result as 

da '\2P b .k s 2p-k s ~ 



The 0(A/E) contributions from Eq. (92) are related to the spin state of M 2 . Because 
up-to twist-3 order, q 2 qz pair for the M 2 can be proportional to ^75, 75, or cr a /3, the 
substitution of these spin terms into the trace over those 0(A/E) terms in the internal 



propagators, such as those last two terms in the right hand side of Eq. (92), results in 
vanishing results by an explicit manipulation. This concludes that the 0(1/ E) soft contri- 
butions vanish and the uncertainties are of 0(1/ E 2 ). In summary, we have shown that the 
V^sih vanishes up-to 0(1/ E) explicitly. The uncertainties are estimated to be of 0(1/ E 2 ) 
which has beyond our precision used in this paper. 
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2. Collinear finiteness 



We now consider the collinear part V^ C {1). It is convenient to combine / M + k c = I'' 1 
TtWwiW M - *^rJ j^yVF* 2^ 2 p ■ k c Vl)] 



x ^[(r^-f^^Cfa)]^. (96) 

We separate the internal parton propagators into the long distance and short distance parts, 
and analyze their contributions, respectively. Let's first consider the term 



Tr[(7°(& + -^)r, - + - J ^)T a )0 M2 «] , (97) 

I 2n-V [I) 2n-l 

where we have explicitly written the long distance term 



^ l(7 F frri W Ta) ^ (1)1, (98) 

and the short distance term 

Tr[( 7 Q -^f,, - t^T)ct> M2 {l)] • (99) 
2n • V 2n-l 

Because the short distance part is suppressed than the leading part of the long distance term 
by an 0(A 2 /E 2 ) factor, we may safely ignore the contributions associated with the short 
distance term. 



In Eq. (98), the 7" can be r y a = fin a , fin a , or, 6^,7"' . The final result also depends on the 
spin state of the M 2 . Because up-to twist-3 order, g 2 <?3 pair for the M 2 can be proportional 
to 0,-75, 75, or e±cr a j3, the substitution of these spin terms into the second trace term Eq. (96) 



leads to the following nine results. We explain them term by term: 



7° = fin a and [9293] ^ 07s : the trace term in Eq. (98) vanishes due to n = 0. 



7° = fin a and [92^3] oc ^75: the trace term in Eq. (98) is proportional to n a . The 



contraction of n a with the first trace term in Eq. (96) gives 



1 jV 2P 6 -fc c 2p-k c J P 3 

This is because n ■ Pb = n ■ p and P& • kc — p ■ kc + 0(A 2 ). The errors are from the iji 
component of fk c . 
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7° = da^ a ' and [52^3] oc $75: the trace term in Eq. (98) 



^Tr[( 7 a '^r j -f l( | 2 7 a >7 5 )] C x^ 



(101) 



is proportional to d°,. The leading contributions of the trace term in Eq. (98) are of 
0(E/A 2 ). The contraction of d°, with the first trace term in Eq. (96) results in 



«x T rp / r ( 2P b,a - lafic) (2Pa+|c7a) F 



2P b • k 



a 



2p-k 



c 



oc 



(102) 



p ■ &c 

which is of order 0(A/E 2 ). Due to the loop integration over fee, the single hex factor 



selects another kc,± factor from the second trace term in Eq. (96). The product of 
these two trace terms are of 0(A 2 /E 3 ). 



7° = fin 01 and [9293] oc 75: the trace term in Eq. (98) is then proportional to 

, 1 1 



'n-l' n .l 



(103) 



The contraction of n a with the first trace term in Eq. (96) gives 



n"Tr[f i (r i g%^^ - 8^±gg^ r< )] 



2P 6 • fc c 
oc — — (P& -n-kc-n) = 



2p-k 



c 



a 

aE 



(104) 



where we have used the fact that Tj £g> E; should be — 2(5 — P) (g> (S + P) and the 
parameterization fc^- = aEn^ + ■ ■ ■. The combination of these two trace terms are of 
0(1/E 2 ). 



7° = fin a and [9293] oc 75: the trace term in Eq. (98) vanishes similar to Eq. (94). 



7 C 



d" 7° and [92^3] oc 75: the trace term in Eq. (98) is proportional to k% ± and is 



of order 0(A/E). The contraction of k^ ± with the first trace term in Eq. (96) gives 



contributions of order 0(A/E 2 ). The product of these two trace terms is of order 
0(A 2 /E 3 ). We neglect it. 



7 C 



fin a and [9293] oc e^a^ u : only Ti ® Tj = —2(5 — P) <g> (S 1 + P) can contribute. 



The second trace term in Eq. (98) gives a result proportional to 

, 1 



+ 



1 



"n-l' n-l' 
34 



(105) 



which is of order 0(A/E). The contraction of n a with the first trace term in Eq. (96) 
gives a contribution of 0(1 /E). The combination of these two trace terms is of order 
0(A/E 2 ). 

7° = fin a and [^2^3] oc e^/a^: only Ti <g> = —2(5* — P) <8> (S + P) can contribute. 



The trace term in Eq. (98) is proportional to n a . The contraction of n a with the first 



trace term in Eq. (96) gives a vanishing result with errors of 0(A 2 /E 3 ). 

In summary, the Tr[Vc(/ 2 )0 M2 (O] vanishes up-to 0(1/ E). Similarly, we can show that the 
term TrfV^cKO^ 2 ^)] vanishes up-to 0(1/ E). Therefore, up-to 0(1/E), we can neglect 
collinear divergences. Note that the above estimated errors should be multiplied by a factor 
A 2 E for the [9253] m 2 s pi n state being proportional to ^75, and A 2 hm 2 f° r the 75 an d • 
spin states, respectively. 

3. Infrared finite one loop vertex corrections 

Based on the above analysis for the soft and collinear limitings of V(l), we now describe 
how to use the subtraction method to show the infrared finiteness of Tr[V(l)(fi M2 (l)]. The 
proof is given by the following series of identities 

dH 



Af = Ff^(q 2 ) J ^Tr^(O0 M2 (O] 

= Ff^(q 2 ) J 7 ^ I (Tr[(V^,(/)0^(/)] + 



F r Ml (q 2 ) / T^imviAi)^)} - Tr[^ ;5 (/v //2 (0] - nvir, c (m M2 {i)\ 



(2tt)< 

dH 

(27) 

-Tr[V ij]C 4)<t> M2 (l)\) + Tr[^, ;5 (/)0 M2 (/)] + Tr[V lJ;C (1)^(1)} 

+Ti[V ll ,c4)<P M2 (l)] + ---] 
Ff^(q 2 ) f 1 dv(Tr[V i:i ., H (v)<t> M *(v)] + •■•), (106) 



3 

•> u 

where the Vi^n(v) is defined as 

V ij;H (v) = Vijiv) - V ij;S (v) - V ij;C (v) - Vw(v) , (107) 
and the dots represent higher order terms which are of order 0(1 /E 2 ). In the last identity, 



the integration transformation 



/ 



d l 4> M2 (l)^ f l dv<f> M *(v) 



(2tt) 4 ^ v ' Jo 
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has been performed. Up-to twist-3 order, the spin factorization can be obtained by using 
the identity 

(108) 

where e± ■ a = <J a (}i and e°f = e a ^ r]1 n r) n 1 . The remaining task is to complete the one 
loop momentum integration over k by using the naive dimensional regularization. The final 
expression will be given in Section |Vj 



H. Collinear expansion for corrections from penguin contractions 



The summation of the contributions from the two diagrams in Fig. [3] (e) and (f ) contain 
the following expression 

Fr'Hq 2 ) J j^yjmm m (l)} (109) 
where the function Pij(l) is defined to contain the quark loop 
Tr[P,,(/)^(/)] 

+ w^^4^ 4tt^-^ • (110) 

The internal quark loop momentum k and the internal quark mass m q has been used in the 
above expression. To extract the leading contributions, we perform the collinear expansion 
for Pij(l) in the following way 

f)P 

Due the momentum derivative over /, the second term in the right hand side of the above 
equation is suppressed by a factor of 0(1/E) than the first term. We first concentrate on the 
first term P^-^). Because both the recoil and emitted mesons are energetic, the radiative 
gluon's momentum scales as (p + I2Y ~ (E,E,0). The infrared structure of Pij(l) then 
depends on the quark loop integration over k, in which k could be soft or hard. If k is 
soft, the numerator of the integrand of the k integration decreases one power of k 2 . With 
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(a) lb) 

FIG. 5: The hard spectator diagrams. 

respect to the loop integration over hard k, the loop integration over soft k is suppressed as 
0(1/ E 2 ). The hard k region is dominated. The integral transformation is made to transform 
(j) M2 (l) into 4> M2 (v). The k integration inside Pij{v) is calculated by using NDR. The explicit 
expression for the above penguin corrections is given in Section |Vj 



I. Collinear expansion for hard spectator corrections 

The relevant expression for the spectator diagrams depicted in Fig. [5] is written as 

ImJwfl ^w B MMi B )^(h)^)] , di2) 

where the spectator scattering function H(l B , / 1; l 2 ) is expressed as 



Tr[H(l B , h, 1 2 )MIb)<P Mi (/i)0 Ma (k)} 
-Tr[r i( M/ B )7 a Ml (/i)] 



(113) 



where k = l B — h- We first perform the collinear expansion of H(l B ,li,l 2 ) with respect to 
the collinear momenta h, % = B, 1, 2, 

H(l B ,h,l 2 ) = H(l B ,h,l 2 )+Y: -EjU-\ k= t l (U-U) v + "- • (H4) 

i=l,2 OL i 

Since the M\ and M 2 mesons are energetic and the spectator quark of the B meson can only 
carry soft momentum, the virtual gluon momentum k can be soft or hard-collinear. For soft 
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k , the second trace term vanishes 



Tl[i ^k Tl ~ f ^)^ M2 ^)] * + °(V^ 2 ) , (115) 

where the errors are estimated to be of 0(1/E 2 ). This shows that the hard spectator diagram 
contributions are free from infrared divergences up-to 0(1/ E). However, it has been noted 
in the literature that there exist end point divergences Xjj at 0(1/ E) as mentioned in 
Introduction. The term X H becomes divergent because the pseudoscalar LCDA ^^(u) is 
a constant. Since the constant model for ^^(u) is determined by the equation of motion 
Eq. Q, as we have shown in Section [IlJ the end point divergences can be identified as a 
mixing effect between the twist-3 and twist-4 LCDAs for the M\ meson. In addition, we 
also showed in Section [n] that, in the energetic meson limit, the pseudoscalar LCDA <f)p 4l (u) 
is reduced to be ^^(u), which is no longer a constant according to the reduced equation 
of motion Eq. (5). The divergence in Xh as u — > is then regularized by the The 
explicit expression for this fact will be given in Section |Vj 

We now show that the 0(1/ E) contributions can only come from the two parton twist-3 
LCDAs for the M\ meson, and the similar contributions from the two parton twist-3 LCDAs 



for the M 2 meson are vanishing at 0(1/E). The 7" in the second trace term in Eq. (113) 
can be 7° = fin a , fLn a , or, 7° . For the spin state of M 2 being 75, the Tj <g> f j can only 
be —2(5" — P)(S + P). For H^(Ib,Ii,12), the dominant contributions in the second trace 
term can only come from tf 2 — tfi, or, ]/ 2 — fi, which are proportional to fi or fl. This selects 



7° = |4ra a , or fin a . For 7" = fin a , the contraction of n a with the first trace term in Eq. (113) 
gives 

n a Tr[r,0 B (Z B ) 7a M ^i)] . (116) 

The contact of ji with <j) Ml (h) extracts one short distance propagator 

— ifi 
2n ■ h 

and a vertex — 27^. The result appears as 

n Q Tr[r,0 B (/ B )7 Q ^(-^X0a /l/3 (/i,/i)] , (117) 

where wt, = gg, — wPnpi. Because the short distance propagator is of 0(1/ E), the dimension 
of hard scattering function is then decreased by one order. The related contributions are 
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FIG. 6: The diagrams for annihilation corrections from final state emission gluons. 
of next twist than what we have considered. The other possibility is that 7" = fin a , which 



which is obviously vanishing. 

For the spin state of M 2 being e± ■ a, the Tj (g) f j can only be —2(S — P)(S + P). Similar 
to the situation for the 75 spin state of M2, the possible twist-3 contribution can only come 
from 7" = fin a , since the fi, v indices of a^ v of e± ■ a are transversal. The second trace term 
is then proportional to 



which is also vanishing. 



J. Collinear expansion for annihilation corrections from final state emission gluons 



then selects the I/2, or f 2 parts of the propagators in the second trace term. The result is 
proportional to 






(119) 



The relevant term for the annihilation diagrams depicted in Fig. [6] , in which the gluons 
are emitted from the final state mesons, is written as, 




(120) 



where 



^[A^l x ,l 2 )^ m {k)4> M2 {h)] 




(fc-fl f f (ft -ft 




1 
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(a) 



FIG. 7: The diagrams for annihilation corrections from initial state emission gluons. 

with k — l\ — l 2 . The collinear expansion of A/(Zi,Z 2 ) is 

A f j (l 1 ,l 2 ) = A f j (l 1 ,i 2 ) + -- - , (122) 

where dots denote the terms of higher than twist-3. Since Mi and M 2 are assumed to move 
in opposite directions and carry energetic momenta, the gluon momentum k can be hard or 
soft. For soft k, the following vanishes 

Tr[(^-F t - f^)0 Ml (/i)7> M2 (^)] ~ + 0(1/ E 2 ) . (123) 
q ■ k p ■ k 

This is because 

Ml (/ 2 )# M2 (/ 2 ) = <p Ah (hW h (l 2 ) ~ + 0(1/ E 2 ) , (124) 
where we have used the property of the long distance part of the parton propagators. 

K. Collinear expansion for annihilation corrections from initial state emission glu- 



ons 



The relevant expression for the annihilation diagrams as depicted in Fig. [7] is given by 



_ ttOsCf %b+1 + m h ) (1/ B + Jt) 

~ ~ 2 ^r Tl[ila (l B + W - rnft " ^ Jh^W la)Mh)] 

Tr[f^ A/ H/i) 7 > M2 (/ 2 )]^ (125) 
with k — h — l 2 . The collinear expansion for A % Mb, h, h) with respect to Ib,i,2 is given by 

A)(l B) hM) = A)(i B MM) + -- - , (126) 
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where dots denote terms of higher than twist-3. The Is is defined as Ib = £,Pb and the £ 
is identified as the momentum fraction carried by the spectator anti-quark of the B meson. 
Some of order 0(A 2 /E 2 ) terms have been added into the propagator of the spectator anti- 
quark in Aj(lB, h, h). This makes the spectator anti-quark be slightly off-shell l 2 B ~ ^ 2 ra 2 oc 
A 2 /E 2 . However, for (V — A)(V ± A) and — 2(5 — P)(S + P) operators, the final expression 
for propagator of the spectator anti-quark appears proportional to l/(u — £) after we have 
cancelled some common factors of the numerator and denominator of the propagator. This 
looks like that the spectator anti-quark carries a collinear momentum ^En^. The subscript 
j in A 1 / mean that j = 1 for (V - A)(V - A) operators, j = 2 for (V - A){V + A) 
operators, and j = 3 for —2(5 — P)(S + P) operators. The superscript i and / in Aj mean 
that gluon emissions start from the initial or final state. The convention used here is that 
the Mi contains an anti-quark from the weak vertex. The anti-quark carries longitudinal 
momentum fraction as u. The M 2 contains a quark from the weak vertex with momentum 
fraction v. Because the weak annihilation terms are power suppressed than the leading hard 
spectator interactions by a relative factor we identify the weak annihilation correction 

as subleading corrections. Within our precisions, only terms from twist-2 DAs for the light 
mesons are considered. Since Mi and M 2 are assumed to move in opposite directions and 
carry energetic momenta, the gluon momentum k can be hard or soft. For soft k, the 
following contraction vanishes 

Tr[( ^l ri ~ r */vi )<M/B)] - + 0{1/E2) (127) 

The hard k contributions should give finite results. However, this requires some cares. 
The second term in the first trace term of li, Iq) is proportional to (u — The 

conventional approach is to neglect the £ as the authors of [17] has done. This then results 
in an end-point divergence Xa as u — > 0. To solve this divergent problem, we propose to 
retain the £ dependence in Aj(ls, h, h)- By introducing a model for </>s(C)> ^ ne divergent 
term Xa appears as 

X A - C d£MZ) t T^F^p 1 W • (128) 
Jo Jo (u — t)u 
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IV. COMPARISONS WITH OTHER EXPANSION SCHEMES 



In this section, a comparison between different calculation schemes for the contributions 
involving the pseudotensor LCDA of the final state light mesons will be given. 

A. BN scheme 

The calculation scheme proposed by Beneke and Neubert in [T7J [2D] is reviewed below 
for comparison. We denote this calculation scheme as the BN scheme. Let H(u,k±, ■ ■ •) 
represent any hard scattering function in the amplitude 

f 1 du ■ ■ ■ Ti[H(u, k±,-- -)® M (u, k±)] , (129) 
Jo 

such as the one loop vertex function V, the one loop penguin function P, the hard spectator 
function H. The twist order of the annihilation contributions involving the twist-3 pseu- 
dotensor LCDA of one of final state light mesons in the B — > MiM 2 decays are identified 
as twist-4, they will be neglected in this comparison. The variable u and k± in H(u, k±, ■ ■ •) 
denote the momentum fraction and transverse momenta carried by the partons, respectively. 
The momentum carried by the quark q is denoted by k q and the momentum carried by the 
anti-quark q is denoted as kq. We write k q and kq as 

k 2 

h£ = uEn» + -^=n" + E 
q 2uE 

k 2 

kg = uEn" + -^-n» - E . (130) 
y 2uE 

The momentum fraction u is defined as« = n-k q /n-P with P being the meson's momentum 
as defined in Eq. ([7]). The u = 1 — u is defined as u = n-kq/n-P. The meson state $ M (u,k ± ) 
is defined as 



■jar- i 



$ (u, k i 



( i f \ ■ -a 3 ( t ) 'ai U ) , ■ v-a^M 

Aip75 I <Pp{u) -ta a/3 n n p — h ia a(3 En — q — ~dk±^ 



(131) 



where <f>p(u), (p^(u), and 4>^f(u) are the twist-2 and twist-3 LCDAs for the pseudoscalar 
meson M, and 0* f/ (u) = d§ a {u)j du. The d/dk±,p in last term of $ M (u, k±) is defined to be 
applied on the H(u, k±, ■ ■ •). After the d/dk±^ has been applied on H( collinear 
limit k± — > is followed to be applied (T71 120]. Because H(u, k±, ■ ■ •) can be expanded in 
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k± , one can use the following transformation [TT] [20 

8 2k? 



dk 



to simplify the calculation 







-k 



(132) 



(133) 



where = g l3X — n /3 n A — n x n^ and the bracket in (k±k±) represents an average over the 



dk ± p 

azimuthal angle in the integration over k±. 

The next step is to use the equations of motion, Eq. Q, to rewrite the twist-3 part of 
$ M (w, k±) in the following expression [T7J [20] 



75^ — -M u ) 



kq ■ k q 



(134) 



tw3 4 

We note that some terms of 0(A/E) have been added to arrive at the above compact form. 



In addition, the solution based on Eq. uAn 4>p{u) = 1 and <j>„{u) = Quu have been used to 



derive Eq. (134). Finally, the Eq. (JT34J) is substituted into Eq. ( |129[ ) to arrive at 
ifpfJ-p rl 



du - ■ ■ Tt[H(u, k±,-- -)75 



rqrq - 
kq ' kg 



4 jo 

For comparison, we suggest to use the following expression 

_ ifp^p 
4 



<j>p(u) . 



(135) 



f du---[ Tr[H(u, ■ ■ -)^\(j) p (u) 
Jo 



o 



+i Tr[ 



dH{u, k q± , 

dkq"± 



-0"a/375j 



(136) 



Using this calculation scheme, the authors of [TTI |2"U] found that the relevant vertex con- 
tributions and penguin contributions are IR finite, while the hard scattering contributions 
contain the divergent term Xh- 

B. DYZ scheme 



Du et. al. j23[ El] proposed to perform the calculations (such as the Dirac algebra and 
the integration over the virtual gluon's momentum) for the contributions involving (f>^(u) 
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in the coordinate space. After the calculations have been completed, the expression is then 
transformed into the momentum space. For comparison, we write the resultant expression 
in the momentum space as 

-^I'du... mH^...) lb } ( p p (u) + ^T,[ ^H ^ (137) 

Using this calculation scheme, Du et. al. found that it is necessary to require a symmetric 
criteria for the <pa{u) to regularize the infrared divergences from the vertex and penguin 
contributions. However, even requiring the symmetric criteria, the IR divergences still exist 
in the hard spectator contributions. 



C. CE scheme 



Based on the analysis given in Section III, the relevant amplitudes for one loop contribu- 
tions with a pseudotensor LCDA calculated under the CE scheme can be formally written 

as 



where 



- P ^ P J du---[ Tr[H{u, ■ ■ -)l^p{u) - ^Tr[H(u, ■ ■ -)e?^]0«r(«) 

1 /"l ' 
+ 8 JO " Tr ^^'"' ' ' ')° r a/375]^0g /3M '(M, u) , 



138) 



JO 2,71 



< "'(M\q{0)a al3 j 5 id^{Xn/E)q(Xn/E)\0) 



(139) 



Under the CE scheme, the analysis given in Section III and the explicit expressions for the 
relevant one loop contributions in the next Section, show that the vertex, penguin and hard 
scattering contributions are IR finite up-to twist-3 under the two particle approximation. 



The last term in Eq. (138) is a of higher twist than three. We retain it in Eq. (138) is only 



for comparisons with the other schemes. 



D. Comparisons 



We now summarize the differences between these three calculation schemes. 
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1. The interpretation of the derivative hard function 

The derivative functions dH(u, l±, ■ ■ -)/dlj_ in the BN scheme or dH(u, I, ■ ■ -)/dl^ in the 
DYZ schemes are interpreted as short distance hard scattering functions related to 4> a {u). 
The derivative d/dl± in the BN scheme, or d/dl 13 in the DYZ scheme come from the co- 
ordinate variable in the spin projector associated with the 4> a {u). In the CE scheme, the 
derivative function dH(u, I, ■ ■ -)/<9/ M arises from a Taylor expansion of the H with respect 
to the collinear momentum I. By a simple manipulation, we can write the following corre- 
sponding relations between the terms calculated in the BN and the CE schemes as 

Tr[H(u, ■ ■ -)a a ^ 5 n a nT N ^ - • • -)ef cr af3 ] CE ^(u) , 

(140) 

Similarly, we can also obtain the following relationships between terms derived by the DYZ 
scheme and the CE scheme as 

n^%r^* a TJi ZE ^fp»P^ - *^[H{u, ■ ■ .)efa aP f E f Pf i P Mu) 

-iTr[H,(u, u, ■ ■ y^^d^f^u, u) . 

(141) 

In the above expressions, we have assumed that (f)' a (u) and (j> a {u) are scheme independent. 
The BN, DYZ, and CE superscripts in the hard scattering functions represent the functions 
having been calculated under the BN scheme, the DYZ scheme, and the CE scheme, re- 
spectively. Under the two particle approximation, H^u.u, • • •) is identical to (dH/dl^)^. 
And, dH(l, ■ ■ -)/(9/ M is 1/E suppressed than H(l, ■ ■ •). The main differences between the CE 
scheme and the other two schemes are from the spin projector ie al3vtJl En a introduced in the 
BN and DYZ schemes. The E factor from the spin projector has the effects to increase 
one order in the dimension of dH(l, ■ ■ -)/dl^. Therefore, the product of dH(l, ■ ■ -)/<9P and 
ie al3utl En a becomes of the same order as of order of H(l, ■ ■ •). On the other hand, in the CE 
scheme, the dH(l, ■ ■ -)/dl^ is not assumed to be associated with a large factor E. 
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2. The twist identification 



The analysis made in Section II indicates that the <p a (u) (or </rJr(-u))) should be identified 
as one twist order higher than the twist order of (j>' a {u) (or 0jj.(ii))) under the energetic light 
meson limit. Namely, 4> a (u) is of 0(A/E) at u = 0(1) and (j>' a {u) = 0(1) at u = 0(1). The 
term 



6 



isl/E suppressed than the term 

-iTr[H(u,- ■ -)a a ^ 5 n- n 



arl 0l(BN) 



6 

Thus, the calculations made by the BN scheme under the energetic meson limit become 
consistent with the results obtained by using the CE scheme. This fact can be seen in the 
next Section, where the explicit results for the one loop contributions calculated by the CE 
scheme will be given. Unlike the BN and CE schemes, the identification of different twist 
order is unclear in the DYZ scheme. The contributions of different order of magnitudes are 
mixing together. 



V. TWIST-3 TWO PARTICLE CONTRIBUTIONS IN B -► ttK DECAYS 

The predictions for the penguin dominated B decay processes under QCDF approach are 
related to the Xh,a terms. We consider the B — > nK decays as an example to illustrate 
how our results obtained in previous sections can be used to improve our understandings for 
these decays. The matrix elements of the effective weak Hamiltonian can be written as, up 
to 0(l/m&), in the convention of [18] 

(nK\H eff \B) £ \^K\T p + T?»\B) (142) 

V ^ p=u,c 

where 

T p = ai(-KK)5 pu (ub) v _ A <g> (su) v _ A + a 2 (TrK)5 pu (sb) v ^ A <g> (uu) v _ A 

+a 3 (nK) £(s6)y_ A ® {qq)v-A + dl(nK) Y^i^v-A ® {sq)v-A 
q q 

+a 5 (nK) ^2(sb) V - A ® {qq)v+A + a 6 (nK) £(-2)(g6) s _ P ® (sq)s+p 
g q 
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+a 7 (nK) J2(sb)v-A ® -e q (qq) v+A 
g 1 

3 

+a 8 (nK) ^2(-2)(qb) S -p ® -e q (sq) s+ p 

q 1 

3 

+a 9 (vrK) J2(sb)v-A ® -e g (gg)y-A 

3 

9 2 



(143) 



where (qiq 2 ) v ±A = gi7 M (l ±7s)<?2 and (gig 2 )<?±p = gi(l ±75)92- The symbol <g> in T p implies 
that the matrix elements of the operators in T p are evaluated according to the factorized 
form (irK\j 1 ® j 2 \B) = {n\3i\B) (K\j 2 \0) or (K\j 1 \B)(ir\j 2 \0). The contributions relate to 
rpann ar j se f rom the weak annihilation contributions with a set of coefficients bi(iiK). 
The expressions for a^K) are written as [TBI ED] 

aftMxMa) = (tf< + °^ X -)N i {M 2 ) 



w±i o^Cf 



4vr z 



v [^(M 2 ) + — F^Ma)] + Pf (M 2 ) , 



(144) 



where the upper or lower signs correspond to the odd % or even i. The superscript p is only 
used for i > 3. The leading order coefficients A^(M 2 ) represent the normalization integral 



of the distribution amplitude 0^f 2 or 



lM 2 



Ni(M 2 



0; 2 = 6, 



1; i = 1-5,7,9,10 



(145) 



The vertex corrections are written as 

-l 



Vi(M 2 



o 



du(f)p 2 (u) 
du<p P l2 {u) 



12 In 



18 + ^(«) 



-12 M— ) + 6 — ^(1 — tt] 
m b 



duff*{u) [-6 + %)]; 



1-4,9,10 



5,7 



6,8 



(146) 



where 



h{v) 



,1 -2m, 
3 I In u — in 



1 - u 



+ 



2Li 2 (u) - In 2 u + 



21nu 



l-u 

2 [Li 2 {u) — (1 + in) In u — (u «-> (1 — w)] 



(3 + 2i7r) In u — (u <->(! — u)) 



(147) 
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where we have employed the naive dimensional regularization (NDR) scheme with an anti- 
commuting 75 for the regularization of the UV or IR divergences arising from the loop 
corrections. The calculations for V^(M 2 ), % = 1 — 5,7,9,10 have been checked by using 
the CE scheme. The results have been found to be consistent with the results derived by 
using the BN scheme. This is reasonable because CE and BN schemes are equivalent at the 
leading twist order. We also note that the Vq^ calculated under the CE scheme are also 
consistent with those obtained by the BN scheme. This is due to the fact that the main 
contributions in V r 6)8 (M 2 ) are from the projection onto the 0^ 2 . The detailed calculations 
for V 6 (M 2 ) under the CE scheme will be given in the Appendix A. The integrations 



f' 1 1 

/ dug(u)(j)p 2 (u) = — - — 3«tt , 
Jo 2 



duh{u)(f)^ 2 {u) = , 



(148) 



lead to the values of the vertex corrections at the scale fx as 



Vi(M 2 ) 



12 In 



-121n( 



37 

Y 
13 



3in 



+ 3i7T 



-6; 



1-4,9,10 

5,7 

6,8 



(149) 



The penguin contributions are given by 

a s C F 



Pl{M 2 ) = 



P?o(M 2 ) 
P£(M 2 ) 



AnN r 



AnN r 



a 



9%N C 

a 
97rN r 



Ci 



4 , rn b 2 
3ln- + --G M2 ( Sp ) 



+ C 3 



8 . m b 
- In — 

3 jJL 



+(C 4 + C 6 ) 



In ^ _ ( n/ _ 2)G M2 (0) - G M2 (s c ) - G M2 (1) 



Ci 



3ln- + --G M2 ( Sp ) 



+(C 4 + C 6 ) 



+ C 3 jjlny + i-G Ma (0)-G Ma (l) 
^ln^ + I - (n, - 2)6 M M - G M2 {s c ) - G M2 (1) 



- 2C|f 
(d + iV c C 2 ) 
(d + iV c C 2 ) 



du(f)p 2 (v 



Gm 2 ( s p) 



4 . m b 
- in — 
3 ji 



4 m 6 2 ~ 
3ln- + --G M2 ( Sp ) 



3 ^/o T 



1 



(M 2| 



- 3C 7 e f 



u 

M 2l 



u 



(150) 
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where rif = 5 denotes the number of light quark flavors. The mass ratios s u = and 
s c = ml/ml are introduced. According to the conventions used in [TB"j [20]. the electroweak 
corrections from C 7 - W are neglected. The functions G M 2 {s)axid Gm 2 ( s ) are defined as 



G M2 (s) 

Gm 2 { s ) 
G{s,u) 



duG(s 
duG(s 



ie, 1 
ie, 1 



—4 / dxxil — x) ln[s 



xfl — x)u] 



(151) 



The effective Wilson coefficients are given by = CV 7 — — C e and C^J = C 8 + C 5 . 



J 8g 



Under the energetic limit, we choose the asymptotic form for both <fi P 2 (u) and (j)^ 2 {u) as 



# a («) 



u 



Quu. The asymptotic models for 0p 2 (w,/i) and 4>^ 2 (u,p) are defined 



as the asymptotic limit \i — > oo for the distribution amplitudes. At finite renormalization 
scale, the distribution amplitudes are expanded into Gegenbauer polynomials 



4>p 2 (u, //) = 6uu 



1 + E«n f2 (^n /2 (2«"1) 



n=l 



(152) 



where the Gegenbauer moments a^f 2 (//) are multiplicatively renormalized. The scale depen- 
dence of the Gegenbauer moments a^f 2 (/i) enters the vertex and penguin corrections at order 
a 2 [18j [20] . In the next-to- leading calculation as we have done in this paper, the Gegenbauer 
moments can be neglected. In this approximation, we then arrive at a further simplification 
Gm 2 { s ) — Gm 2 {s). This results in the following identities 



P|(M 2 ) 
Pi{M 2 ) 



p!{m 2 ) + ^a e// 



In 

P?o(M 2 ) + ^G 



8s 
eff 



(153) 
(154) 
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For practical applications, the Gm 2 (s) are evaluated under the previously mentioned 
approximations as 



GmM 
Gm 2 (1) 
Gm 2 {s c ) 



5 2iir 
3 + T~ ' 

85 r- 4vT 

5 2, 32 

- - - lns c + —s c + 16s 

2 



(155) 
(156) 



3 

+I2s 2 c 



y/1 - 4s c 1 + 2s c + 24s 2 c (2arctanh v / l - 4s 



in) 



1 - 



(2arctanh\/l — 4s c — in) 2 + 



(157) 
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The complete expressions for Gm 2 (0), GW 2 (1) and Gm 2 (s c ) with Gegenbauer moments are 
referred to [IB| l2Dj. 

The scale for the vertex and penguin corrections refers to the parton off-shellness in the 
loop diagrams. The typical setting of the scale is chosen to be fi ~ my,, which are the 
scales substituted in the Wilson coefficients Cj and in the hard scattering kernel T 1 . The 
combination of the scale dependences and the renormalization scheme-dependent constants 
in Cj, Vi(M 2 ) and Pf(M 2 ) give renormalization-group invariant results. 

The hard spectator corrections are given by 



H^M.M,) = B t M ^m B f L du C dv 

AmiMo jo t Jo Jo 



NhM 2 JO £ 



for i = 1 - 4, 9, 10, and 



H i (M 1 M 2 ) = B ^m B fdt^P- fdu ('dv 
AmxM 2 jo £ jo Jo 



u(v — £) x u(v — £) 

(158) 



0P 1 H0P 2 (^) +r M 1 < 1 («)^ 2 (^ 



u(v - f ) x u(v - i) 

(159) 



for i = 5,7, and Hq^(MiM2) = 0. Different from the BN scheme, we have introduced the 
distribution amplitude for the B meson and retained the £ dependence in the spectator 
propagator. Since the (^^(u) = 6uu is no longer a constant, there are no end-point diver- 
gences as in those results derived in the BN [T8| 120] or DYZ schemes [23l |2"4"] . Due to the 
soft scale associated with the spectator quark, the introduction of the B meson's DA 4>b{0 
may need to consider the effects of double logarithms from a overlap of soft and collinear 
divergences associated with the spectator quark. This may need to introduce a Sudakov 
form factor to account for the effects of double logarithms. A two loop analysis for the hard 
spectator made by Beneke and Yang [46J indicates that there exists no double logarithms 
up-to 0(a s ). Therefore, at next-to-leading order, we can neglect the Sudakov form factor 
completely. By using the model for 4>b(0 and the fact that 0p fl (w) = 0p fl (-u), we arrive at 
a simple expression for Hi(MiM 2 ) , which is useful for numeric calculations, 

E t {M x M 2 ) = B ^m B {\ + r*)N B , (160) 
Am x m 2 

where Nh = 18/Ab. The model for 0b(£) is assumed to have the form 

N B ee 

[C 2 + est] 



b(o = r „ :iv- 12 (lei) 
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with Nb = 0.133, €b = 0.005 for Xb = 350 MeV. The Nb and e# are determined according 
to the conditions 



Jo 



o 4 X B 

where the errors are controlled within 1%. This is consistent with the conventional assump- 
tion \b < 600 MeV under the condition 3\b < 4A with A = m^ — m^,. The scale dependence 
in r™ l (n h ) is chosen as fj, h = y/A^jj with A h = 0.5 GeV pH 120] . 

The annihilation corrections are expressed in terms of the following Aj, j = l,---,3 
functions 

-1 r l 



A\ = 7ra s I d&eiO I du<p A P h {u) / dvcpf 1 ' 
Jo Jo Jo 



A' 2 = j ' (iConiQ <luo]'. '{u) <lco)l> 



1 1 

(u — £)(£ — v))v (u — £)uv 
1 1 



(l-(u-f)(f-v))it (v-£)uv 



A\ = A{ = A{ = A{ = . (163) 

The subscript j in A % ^ mean that j = 1 for (V— A)(V— A) operators, j = 2 for (V— A)(V+A) 
operators, and j = 3 for —2(5* — P)(S + P) operators. The superscript i and / in Aj mean 
that gluon emissions start from the initial or final state. The convention used here is that 
the Mi contains an anti-quark from the weak vertex. The anti-quark carries longitudinal 
momentum fraction as u. The M 2 contains a quark from the weak vertex with momentum 
fraction v. Because the weak annihilation terms are power suppressed than the leading hard 
spectator interactions by a relative factor l/m&, we identify the weak annihilation correction 
as subleading corrections. Within our precisions, only terms from twist-2 DAs for the light 
mesons are considered. The non-singlet annihilation coefficients are given by 



h = %CiA\ , 

c 


h = %C 2 A[ , 

c 




^ = %C,A[ , 

c 


bl^^[C 4 A{ + C 6 Ai], 

c 




°3,EW — ^72°9^l 7 

c 


Cf 

K,ew = [CioA\ + C 8 A 2 ] . 

c 


(164) 



We note that the coefficients (pi, b 2 ) corresponds to the current-current annihilation, (63, 64) 
corresponds to the penguin annihilation, and (b^EWi ^a,ew) to the electroweak penguin an- 
nihilation. The end-point divergences as the spectator quark approaching its on-shell are 
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regularized by the momentum fraction £, which is defined as Ib = £-P&- Such a definition for 
£ would make the spectator quark slightly off-shell since l 2 B = £, 2 P 2 ~ \ 2 B , which is of order 
Aq CD . The errors from this off-shell- ness for the amplitude is of order 0(A 2 /E 2 ), which is 
of twist-4. Within our precision, the errors can be neglected. 

Using the same approximations for the hard spectator corrections, we arrive at the fol- 
lowing simplified expressions for the annihilation corrections 

A\ = A\ ~ 18™, {^^- + N A ^J , (165) 

where = —12.1 — 1.7i is related to the model for <Pb{0- The convolution integrations 
appear to be overlapped integrals over £, u and v with 4 >1 p 1 { u ) an d (f)p 2 (v). The 

traditional end point divergences as u — > 1 or v — *■ 1 are then regularized by the B meson 
distribution amplitude 4>b{0- Although the averaged value of the £ is much smaller than 
one, however, we argue that it still can not be neglected. The scale dependence in a s (fi) is 



chosen as jj, h = y/AfJl with A h = 0.5 GeV. 

For penguin dominant B — > 7rK decays, the relevant decay amplitudes under QCD 
factorization are parametrized as the following 



A(B~ -> n-K°) = X p 



K- 2 a ?o) + r x ( a 6~ \ a l) 

EW\ 



A 



ttK 



+ (\ u b 2 + (\u + \c)(h + b* w ))B TK , 

-V2A(B~ -> n°K-) = [X u ax + A p « + a p w ) + \ p r«(a p 6 + a p 8 )]A nK 

3 

+ [\ u a 2 + (A u + A c )-(-a 7 + a 9 )]A K7r 
+ [\ub2 + (\u + \c)(h + bi w )}B nK , 
-A(B° -> 7T+^-) = [A m0i + A P K + a? ) + A p rf K + aS)]4r* 

+[(A„ + A c )(6 3 -i&f w )]fW 
V2A{B° -> vri? ) = -> 7T-/?°) + >/2A(B- -> vr ^) - A(fi° -> tt+JT) (166) 

where \ p = V P bV* s , a, = a^trK), and A p af = A u a" + A c a^. The CP conjugation of decay 
amplitudes are obtained by replacing X p — > X* for the above amplitudes. The A n x, Axn, 
and B n x are defined as 

A,k = t^(m 2 B -ml)F B ^(m 2 K )f K , 
Ak* = i^(m 2 B -m 2 K )F B - K (ml)U, 
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G F 

BnK = I^/b/ic/w • (167) 

For numerical calculations, we will use the following input parameters 

Ag } s = 0.225GeV , m b (m b ) = 4.2GeV m c (m b ) = 1.3GeV , m s (2GeV) = 0.090GeV , 

1^1 = 0.41, 1^/^1 = 0.09, 7 = 70°, t(B') = 1.67(ps) , . . 

(168) 

r(B d ) = 1.54(ps) , / w = 131MeV , = 160MeV , f B = 200MeV , 

F*-+* = 0.28 , F^ K = 0.34 . 

For X u and A c , we take the following convention for their parametrization 

^ = tan 2 6 c R b e-^ (169) 

where 

A 2 

tan 2 9 C 
Rh 



1 - A 2 
A 'V^ 



1-A 2 /2,K 



A = |K.| . (170) 

The value of A is taken as 0.22. 

The values of the NLO Wilson coefficients Cj, i = 1, • ■ • , 10 and the LO C%t and C^J , 
at the scale m b = 4.2GeV ( y/A h m b = 1.45GeV) are given by [18] 

d = 1.081(1.190) , C 2 = -0.191(-0.373) C 3 = 0.014(0.027) , 

C 4 = -0.036(-0.062) , C 5 = 0.009(0.012) C 6 = -0.042(-0.086) , 
C 7 /a = -0.001(-0.013) , C 8 /a = 0.060(0.133) C 9 /a = -1.254(-1.380) , 
do/a = 0.223(0.432) , C^ f ' LO = -0.318 , C 8 e g //,LO = -0.151 , 



(171) 



where the parameters A.3 S = 0.225GeV, m t (m t ) = 167GeV, m b (m b ) = 4.2GeV, M\y = 
80.4GeV, a = 1/129, sm{6 w ) = 0.23, a s {M z ) = 0.118 have been used. By using the 
input parameters and the Wilson coefficients, we list the values of ai, % = 1, • • • , 10, and bj, 
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j = 1, • • • , 3, calculated at the scale m& = 4.2GeV as below 

ai = 0.938 + O.OUi , o 2 = 0.351 - 0.081* , o 3 = -0.011 + 0.003i , 

a\ = -0.020 - 0.015i a\ = -0.025 - 0.006i , a 5 = 0.017 - 0.003i , 

r K a u = _o.048 - 0.016i , r^a% = -0.054 - 0.006i , a 7 /a = 0.060 + 0.004i , 

T%al/ot = 0.077 -O.OUi , r^ag/a = 0.072 - 0.07i , a 9 /a = -1.149 - 0.017i , 

a\Ja = -0.378 + 0.082i , a c l0 /a = -0.382 + 0.088i , r A h x = -0.154 - 0.022i , 

r A b 2 = 0.048 + 0.007i , r A b 3 = -0.003 - 0.0005i , r A 6 4 = 0.019 + 0.003i , 
r A bf w /a = 0.179 + 0.026i , r A bf w /a = -0.073 - O.OUi , 



(172) 



in which 



~ m 2 B F B —(0) 



r A =~^ = - 2 ilZ^ * 0.004 , (173) 



and 



2m 2 

r x = i 1 - n - 1-18(0.76) , (174) 
A m b {m q + m s ) 

where = 1.18 is calculated at the scale m h = 4.2 GeV, and = 0.76 is calculated 

A, A 

at the scale y/Ahrrib = 1.45 GeV. Different contributions to the coefficients a, are given 
in Table I for reference. The hard scattering contributions are dominant for 03,05,07,010- 
The penguin contributions are minor for all Oj, i = 1, • • • , 10. The vertex contributions are 
important except of a^^. 

The branching ratio for a B — > nK decay is given by the expression 

Br(B -» ttK) = Yjr^ \A(B^irK)\ 2 , (175) 

The predictions in the CE column of Table II are calculated according to the results derived 
in this paper. The predictions in BN column of Table II are quoted from the paper [20J. 
The predictions in DYZ column of Table II are quoted from the paper j2I] . We observe that 
the predictions in the BN column are about 1.7 — 1.9 times larger than those in the CE and 
DYZ columns. To understand this difference, we employed the theoretic input parameters 
given in the Table 2 and the Oj and b{ values given in the Table 3-5 of [18] to recalculate 
the predictions for the branching ratios for the B — > 7r.fr decays. The predictions are given 
in the BN' column of Table II. To check the consistency between our calculations made in 
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TABLE I: Different contributions to the coefficients m, i = 1, • • • , 10. 



N c 



a s CFCi±i 
47r7V c 



Vi 



Tra s CFCj±i 
AT2 



H; 



P 



total 



ax 



1.017 0.028 ±0.014?- 



-0.107 



0.938 ±0.014?- 



0.169 -0.158 -0.081s 0.340 0.351 - 0.08H 



A3 


0.002 


0.005 + 0.003?' 


-0.018 





-0.011 + 0.003i 


a\ 


-0.029 


-0.002 - 0.00H 


0.008 


0.003 - 0.014s 


-0.020 -0.015? 


a% 


-0.029 


-0.002 - 0.00H 


0.008 


-0.002 - 0.005? 


-0.025 - 0.006i 


a 5 


-0.005 


-0.002 - 0.003?- 


0.024 





0.017- 0.003? 


"6 


-0.039 


-4.27 x 10~ 4 





-0.002 - 0.014? 


-0.041 -0.014? 


a§ 


-0.039 


-4.27 x 10" 4 





-0.007 - 0.005? 


-0.046 - 0.005i 


a-j/a 


0.019 


0.003 + 0.004? 


0.038 





0.060 + 0.004? 


ag/a 


0.060 


-4.75 x 10~ 5 





0.005 - 0.012?- 


0.065 - 0.012?' 


ag/a 


0.060 


-4.75 x 10~ 5 





0.001 - 0.006? 


0.061 - 0.006? 


ag /a 


-1.180 


-0.033 - 0.017?- 


0.064 





-1.149-0.017? 


aio/a 


-0.195 


0.184 + 0.094? 


-0.395 


0.028 - 0.012?- 


-0.378 + 0.082?- 


af /a 


-0.195 


0.184 + 0.094? 


-0.395 


0.024 - 0.006? 


-0.382 ±0.088?' 



TABLE II: Predictions for CP-averaged branching ratios (in unites of 10 6 ) for B — ► TtK decays. 





mode 


CE 


BN' 


BN[20] 


DYZ[21 


datap3] 


B~ 




10.1 


12.1 


19.3 


12.6 


23.1 ± 1.0 


B~ 


-» TT°K- 


6.4 


5.3 


11.1 


6.7 


12.8 ±0.6 


B° 


-> TT+K- 


9.6 


10.2 


16.7 


9.1 


19.4 ±0.6 


B° 




3.9 


6.0 


7.0 


4.3 


10.0 ±0.6 



the BN' column and those in [18], let's make an example for comparison. For example, the 
prediction for Br(B~ — > -k~K°) decay is calculated to be 12.1 x 10~ 6 , which is close to the 
central value given in the Eq. (93) in the paper [18] 

10" x Br{B- - n-K°) = [14.1^; 4 (m s )^(X A )] x 



F B- 



-(0) 



0.28 



(176) 
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where the first error is due to the parameter variations shown in the Table 2 in that paper, 
and the second error is from the uncertainty due to power corrections from weak annihilation 
and twist-3 hard spectator contributions. We find that our calculation is consistent with 
that made in |18j . We argue that the difference between the predictions given in the BN' 
and BN columns maybe due to different methods used for determining the central values of 
the predictions among many sources of theoretical uncertainties. In average, the predictions 
made by the CE, BN' and DYZ schemes are only about one half of the experimental data. 
The predictions made by the BN scheme given in the paper [20] are consistent with but 
lower than the experimental data [71]. 

VI. DISCUSSIONS AND CONCLUSIONS 

In this paper, we have shown that the factorizability of the QCDF amplitudes for B — > 
PP decays with twist-3 two parton contributions can be preserved under the energetic 
meson limit that the two final state light mesons carry energetic momenta. Namely, we 
have extended the factorization theorem Eq. ([I]) to 0(a s ) and 0(l/mb) under the two 
parton approximation. The factorizability is shown by the following facts: (l)Under the 
energetic meson limit, the pseudoscalar distribution amplitude for a light pseudoscalar meson 
is allowed to be non-constant by the equations of motion for the quark. (2)The non-constant 
4> p (x) is then used to regularize the end-point divergences in the hard spectator corrections 
at twist-3 order. (3)By retaining the dependence of the momentum fraction variable of 
the spectator quark of the B meson, the end-point divergent problem for the annihilation 
corrections at twist-3 order are solved. Based on the factorization for the matrix element at 
the twist-3 order, we have constructed a collinear expansion (CE) scheme for calculations 
of the hard scattering kernels of order 0(a s ) and 0(AQ C D/ m b)- The results were applied 
to make predictions for the branching ratios of B — > nK decays. The predicted averaged 
branching ratios of B — > irK decays are only about one half of the experimental data. 
Because the end point divergences in the hard spectator and annihilation corrections as found 
in previous studies have been regularized, the strong phase related to these two contributions 
is predicted to be universal as 4>a = 8.2°. 

The predictions for the averaged branching ratios of B — > 7iK decays made by three 
schemes, the CE, BN' and DYZ schemes, all have similar magnitudes. Although the CE 
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scheme contains no significant improvements in phenomenology than the other two schemes, 
however, it has reduced large uncertainties in the Xh,a terms, theoretically. In literature, 
the Xh,a terms are modeled as 

X HA = \n T ^ B (l + p HA e^), (177) 

where Ph,a < 1 an d 4>h,a are unknown parameters and process dependent in general. In 
order to reduce the theoretical uncertainties , pu,A and 4>h,a are further assumed universal 
and their values are determined phenomenologically [20]. This reduces the prediction power 
of QCDF. On the other hand, in our approach, the relevant terms are finite and their values 
are calculable. The introduced B meson DA <Pb(0 in the weak annihilation and twist-3 
hard spectator contributions results in no additional uncertainties than the other schemes, 



because the parameters in the </>_b(£) Eq. (161) are completely determined by the first two 



moments of the 4>b(0 Eq. (162). Since the most contributions of A\ 2 are from the end-point 
region, it is also interesting to check whether different models for the 4>b(Q can result in 
different values for A\ 2 . As found in [75], the end-point behaviors of the model for the 



4>b(0 are well controlled by the first two moments Eq. (162) and almost independent of the 



parameterization form of the assumed model. This implies that the A\ 2 given by Eq. ( 161 ) 
are almost model independent. In summary, the predictions under our approach are not 
only free from the end-point divergences but also independent of the model for 0# (£). In 
addition, we have unproven the theoretical uncertainties to be of order 0(1 /ml). 

We emphasize that the methods proposed in this work for resolving the end-point diver- 
gences is not ad hoc but general. For example, in the charmless B decays with one scalar 
meson in the final state, the chirally enhanced corrections are necessary. The twist-3 two 
parton DAs for the scalar meson can get involved. Similar to the pseudoscalar meson, the 
equation of motion for the twist-3 two parton DAs for the scalar meson are also used to 
determine the DAs [76] 177] . The physical situation in these decay processes is similar to 
that in B — > nK decays. Therefore, the analysis given in Section |TT] and the calculation 



scheme given in Section III can be used. We will show this fact in another place. Another 
important contribution of this paper is that the proposed CE scheme not only provides a 
systematic method for including higher twist contributions, but also is consistent with the 
QCD factorization scheme. This is the first method in the literature that can systematically 
calculate the higher twist contributions within the factorization approach. Last, the com- 
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plete twist-3 contributions need to consider the three particle DA. We plane to study this 
issue in another preparing work. 
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APPENDIX A: CALCULATIONS OF TWIST-3 VERTEX CORRECTIONS 

In this Appendix, the calculations for the one loop vertex contributions to the matrix 
element of operators Qi, % = 6 will be present. The amplitude for the Feynman diagram as 
depicted in Fig. |3|a) is written as 

/ r»(S)(y-A)(y+A)\vertex,l-loop 
W5 /(a) 

= -^vr~ y rfMr(3) I dx Jo dy ij^r [k^2( X p b+yuq ).kr 

x (M 1 \q-f(l + j 5 )M™ 2 {u) - l -e ± ■ *^{u))r{ui - m»(l - 75) (A - * + m 6 ) T «6|S> . 

(Al) 

where NDR has been used. In the above expression, we have used the following spin state 
expansion for the matrix element (M2\q(0)q(Xn/ E)\0) 



infty rj\ if** / I ~ \ 

-e-^(M 2 \q(0)<l^n/E)\0) = -^[ 7 5#^(«) + A*? a - ^ " ^»(«)J] 



(A2) 



Perform a substitution of — > + xP b + ywg for k to arrive at 

/ n ( v ~ A ) ( v+ A ) \ vert ex > 1 ~ 1 °°p 

\V5 /(a) 



ifhhUM ™* ^ J q du(f)p h {u)T(3) dx dy J 



(2n) D [k 2 - (xP b + yuq) 2 } 3 

1 

2 

(A3) 



x(M 1 |g 7 ^(l + 75 ) ( 75 - - ei • a) 7 a (^ - xft - #) 7 „(1 - -f 5 )(xP b - yuj - $ + m b ) 7a 6|5) , 
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where we have used <fip l2 (u) = 4>^ 2 (u). The contributions with odd number of k give van- 
ishing results. Completing the loop integration over k by NDR gives 

/ / ^ l (8)(V-A)(V+A)\vertex,l-loop 
W5 /(») 

na s C F ^ D /-i , M ri f-* (-1) 3 T(3 - D/2) ( 1 \ 3 ^ /2 

= «/M 2 /iM 2 77 / 2 (w)r(3) dx dy 



N c Jo ^ p v 7 v 'Jo Jo tf (47r)^/ 2 T(3) \(xP b + yuq) 2 , 
1 
2 

x- 1 1 

(Mil^^l + 75 )( 75 - -e ± ■ a) 7 a 7r;7/ ,(l - 7 5 )7A7a6|S) } . 



x { {M x \q^{\ + 75 ) ( 75 - 2 e ± ■ a)i a {uyi - xf b )^{l - 75) {xf\ - yui + m b ) 7Q 6|5> 
g vX T(2 -D/2) ( 1 



2 r(3 -D/2) \(xP b + yuq) 2 



(A4) 



We note that the first term in the bracket of Eq. (A4) is finite under x, y — > 1 while becomes 
divergent as x, y — > 0, and the second term in the bracket of Eq. (A4) is finite under x, y — > 
while becomes divergent as x, y — *■ 1. This implies that the first term is infrared divergent 
while the second term is ultra- violate divergent. We apply D = 4 + 2a for the first term and 
D = 4 — 2e for the second term. 

The Dirac identities in D-dimension 

7% = D , 

rri, = —{d - 2)y , 

yvy 7/l = 4^ - (4 - d)t^ , 

7 M 7 a 7 /3 7 7^ = _2 7 7 7 « 7 /3 + (4 _ j D) 7 « 7 /3 7 7 ; ( A5 ) 



can be used to simplify the Dirac matrix. After completing the spin algebra, we found that 
only 7 5 0^ 2 (m) can contribute. The UV part of (Q^ V A ^ v+A ^y^ tex ^~ lo °P j s 



/ n (8)(V-A)(V+A)\vertex,l-loop 
\V5 )(a),UV 

-2if M ^(M x \q{\ - l5 )b\B)^-£ j~ du$*(u) 

x(l - -e) ( Tie) ^ dx dy— — — . (A6) 

v 2 J \ m 2 J K ' Jo Jo y [x(x + yu)Y V ; 



The integrals over x and y are 



/ , 2 uhiu 

dx / tZy-p-p — = 1 + 3e H — . (A7) 

Jo [x(x + yu)\ e u 
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By substituting the above identity, the UV contribution becomes 



/ n (8)(V-A)(V+A)\vertex,l-loop 
W5 )(a),UV 



= -2if M ^ M2 (MMi-^)b\B)^£j o du^\u) 

,1 , u 2 1 «ln«. . , _ 

x --7 £ + 4vr + ln^2 + - + ^^ . A8 
e mi 2 u 

The calculation for the IR part of (Q^)(v- A )(v+ A )y^rtex,i-ioo P re q U j res some cares. As we have 
shown in Sec. Ill, the different components of the vertexes 7 M , 7°, 7^ lead to contributions 
of different magnitudes. The IR part of (Q(v- A )( v + A )y£rtex,i-i oop appearg ag 

/ n (8)(V-A)(V+A)\vertex,l~loop 
W5 /(a) 



= if M2 m™ sC ^ f du^(u)r(3) £ dx f Q X dyj^ 



(2ir) D [k 2 - 2(xP b + yuq) ■ k} 3 
x (M 1 |^(l + 75)75(2(yM - - x-ffi^l - -f 5 )(2xP b a + (xm b - yui) la )b\B) . 

(A9) 

By some manipulations, the expression becomes 

/ n (V-A)(V+A)\vertex,l-loop 
W5 /(a),IR 

1 

= 2i/ MaAiAf2 (M 1 |g(l- 75 )6|S)^^- jf rf M 0f 2 («) 

x (^f] r(l - a) rfx dy . . NlR ... , (AIO) 

where 

N IR = [2Duxy + 4uxy - (4 + 2D)xx] . 
After completing the integrations over x and y, the ^( 8 )( y ~ A )( y+A )^ertex,i-/oo P ^g^^gg 

/ n (V-A)(V+A)\vertex,l-loop 
W5 /(o),/fl 

= 2^ M2/ i M2 (M 1 |glf-75)&|5)^(^ 4 ^ 2 ^ r(l-a) £ duff*(u) 



r l 21nw 1 , 9 . 1. _ 31nM n /,„,s 

x — + + ln-u + In -u — 2L22(1 ) + 1 H z — . All 

a 2 a 2a u u 

As a result, (q^-^)^))^,!-^ ig written &g 



/ / ^ ) (V-A)(V+ J 4)\Dertex,l-ioop 
\V5 /(a) 
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-2if M2 n M MiW - ^) b \B)^f I' du$\u) { (- - lE + 4vr + In 4 + \ + ^) 

AttJM c Jo e m£ 2 u 

/47r/i 2 \ a . , r l 2\nu 1 , , 9 , 1, 31nu ni 

- — V T l-a - + + — -lnw + ln 2 u-2Lz 2 (1 - - +1 + } . 

\ m( J a 2 a 2a u u 

(A12) 



The calculations for the Feynman diagrams depicted in Fig. |3Jb-d) can be done in a 
similar way. The amplitudes for these three diagrams are written as 

i n {V-A)(V+A)\ vertex, 1-loop 
W5 /(&) 

= 2zf M2f x M2 (M^? l-7 5 & B /-f / ci^ /2 ( M ){(--7i=; + 47r + ln^ + - + ) 

47riV c jo e m£ z u 

{ Airu 2 \ a _ ,„ . r 1 2 In u 1 2 . 1 3 In -[L .. 

- — ^ T 1-a — + + — -lnM + ln 2 M-2L« 2 1-3 + ! + } • 

\ mi J a 2 a 2a u u 

(A13) 



/ /O A) (y +^) \ vertex, 1- ioop 
\V5 /(c) 



2</ WaA i A r !1 <M 1 |g(l - /' du^(u) { ( 1 - lE + 4vr + In 4 + £ - ton + iir) 

A-kJSc Jo f e ml 2 

/47r/i 2 \~° w , u 2 21nw 3 + 2m „ n . 4 . l2 27 + 27T,, 

- — ^~ r i-a — + 4 + 2«7T lnu + 4zvr + ln 2 M + } . 

\ mi J a z a 2a 6 

(A14) 



/ D ( V ~ A ) ( V + A ) \ vertex, 1-loop 
\V5 1(d) 

-2if M2 mMiW ~ ^)b\B)^- f 1 du^ h (u) { (- - lE + 4vr + In 4 + \ ~ h« + in) 

Airi\ c Jo e mf 2 

(A-Ku 2 \~ a . Nr 2 21nu 3 + 2m , 2 27 + 271%, 

(A15) 



The summation of the above contributions from Fig. [3] (a)- (d) gives the in Eq. (88). It is 
obvious that the IR divergences are cancelled. The UV divergences are regularized by means 
of the dimensional regularization and the remaining scheme constants are subtracted by the 



MS subtraction scheme. Our calculation is consistent with the result derived by using the 
BN scheme [20J. We also note that, for the vertex contributions, the pseudotensor LCDA 
does not contribute in both CE and BN schemes. In the DYZ scheme, the pseudotensor 
LCDA can involve in the vertex contributions and there are associated IR divergences . 
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The authors in [23] argued that a symmetric parameterization for pseudotensor LCDA is 
necessary for eliminating the associated IR divergences. Although both CE and BN schemes 
obtain similar results for the vertex contribution Vq(M 2 ), the meanings for the pseudoscalar 
LCDA ((f)p l2 (u) in the CE scheme and the (j)^ 2 (u) in the BN scheme) are different. In the 
CE scheme, the 0^ /2 (u) = 6mm is used. On the other hand, in the BN scheme, <p^ 2 {u) = 1 
was used. 



APPENDIX B: CALCULATIONS OF TWIST-3 PENGUIN CORRECTIONS 

The calculations for twist-3 penguin contributions from the penguin contractions as de- 
picted in the diagrams in Fig. [3] (e) and (f) are given in this Appendix. We denote the 
penguin topology in Fig. [3] (e) as type-/ and penguin topology in Fig. [3] (f ) as type-//. The 
contributions from Qi 3 appear in the type-/ penguin topology and the contributions from 
(^4,6 are identified as the type-// penguin topology. For the type-/ penguin topology, the 
amplitude from operators Qi,3 are written as 

//~)(V—A)(V—A) \ penguin, \— loop 
Wl,3 / 

C F 7ia sf i A ~ D £ Mo /-i /-i r d D k 1 



N c Jm '^ x Jo Jo J (2tt) d [k 2 - R 2 } 



2 



x ( p - uq y ^ M ^ a ^ 2 ^ ~ \*± ■ <# a («)b M (i - 75) 

(f, - xty - u$) + m q ) la {lfc + x(i> - ui) + m q )j„(l - i h )b\B) , (Bl) 

where R 2 = m 2 (s q — xxu) with s q = m q /m|. For the type-// penguin topology, the relevant 
amplitude from operators Q^q is expressed as 

I r\{V~ A)(V— A) \ penguin,!— loop 
W4,6 / 

C F *a.f- D , ,. M2 ri ri _ f d D k 1 



N, 



c 



/M 2 uf I duJ Q dx J {27r)D _ R2] 



2 



1 .... XM 2 f„ l \ _ 1 . „2M 2 , 



x (p _^ )2 (M 1 |g 7 Q (75< 2 H - 2 e± • ^(^"(l - l,)b\B) 
x 2 { Tr[ 7Q (^ - xty - ui) + m q )^(l - 75 )(# + x{i> - ui) + m q )\ } . (B2) 



q=u,c,t 



where R 2 = m 2 (s q — xxu) with s q = m 2 Jm 2 . We first complete the k integration to obtain 

/ /~\{V — A){V — A) \ penguin,!— loop 
\Vl,3 / 
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tfM ^ 2 ^kIo du Io dx 



1 rl 



x { ^R 2 ({Mi\ g7 a (75 (j)™ 2 (u) - 2 e ± • a^{u))^{l - 75)7^7,(1 - 1s)b\B) 
+XX {^nf) J^( M ^ a ^ 2 (u) - \e ± ■ a^(u)mi - 75 )((^ - H) + ™ q ) 



- u{) + m 9 ) 7M (l - 7 5 )6|5> } (B3) 



1 

(p — wg) 2 



and 



\ penguin,!— loop 
\V4,6 / 

ifM ^^jl du ii dx 

x { ^^ (Mi| - 7 a (7 ^M 2(M) _ l c± . a ^M 2{u)mi _ 75 ) 6 | jB)Tr[7a7 A 7M(1 _ 75 ) 7 ,] 



+XX i^mf) ^ M ^ a ^p 2 («) - 2 £ ' ^M)^ 1 - 



x( 2 1%^ - ttfl + mgh/X 1 ~ 7s)((?* ~ ui) + m g )] ) } , * .(B4) 



<j=m,c,6 

where .R = ^fs q ~ zr ^cxu. After completing the spin algebra, we arrive at 



(p — uq) 2 



jo 



dxxx 



[ 7s + ln(4vr) + In - 2 - 2 In Js q - xxu] (B5) 

C TYlu 



and 



1 u 2 
^ [-- lE + ln(4vr) + hi ^ - 21n Js, — xxu] . (B6) 

q=u,c,b o 

After completing the integral over x, we obtain the contributions from different operators 
contained in the Pq(M2) function. Similar to the vertex contributions, the pseudotensor 
LCDA does not get involved in the penguin contributions. The penguin contributions cal- 
culated by the CE scheme are consistent with the results calculated by using the BN scheme 
|2Uj . while are different from the results calculated by the DYZ scheme [23J. 

The calculation for the contributions from the chromo-magnetic dipole operator Q 8 g is 
straightforward in the CE scheme. The amplitude is written as 



i Z\r ^(M 1 |g 7 a (75CH-o e ^C 2 («))[^-7^](l + 7 5 )&| J B)- 



(B7) 
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After completing the spin algebra, we arrive at 



J M fT S tdu(MMUT («) + V€Ku)){l - lB )b\B) 



2% N c . Jo -'^2" p v ' 2 

a s C F r 1 



= (-4if M2 ^ M2 )^£ J o du^{u){M x \q{l - lB )b\B) , (B8) 

where the second line is obtained by using the equations of motion Eq. ([5]). For comparison, 
we note that the calculation for the contributions from the chromo-magnetic dipole operator 
Qsg by the BN scheme gives the following expression [20] 

= {-4if Ma fi Ma )T%- t ducj> p {u){M x \q^-l^\B) , (B9) 



4<kN c Jo 

where the second line is derived by using the Eq. Q. We note that the results from the 
penguin contractions calculated by the DYZ scheme are quite different from our results. The 
detailed expressions for the penguin contributions by the DYZ scheme are referred to their 
original paper [23]. 

APPENDIX C: WILSON COEFFICIENTS 

The NLO Wilson coefficients used in this paper are summarized below for reference. The 
solution to the renormalization group equation for the Wilson coefficients C\, ■ ■ ■ , C±o can 
be written as 



47T 47T 47T \a s {j2) 



C{M W ) (CI) 



where the U$ and Rq are the LO parts and Jo and R\ are the NLO parts of Cj, respectively. 
The a s is used its NLO expression. The expressions for the Uq, Jo, -Ro and R\ matrices can 
be found in |73|. 
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